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1 Half a century of large-scale magnetic fields 
1.1 A premise 

The content of the present contribution is devoted to large-scale magnetic 
fields whose origin, evolution and implications constitute today a rather in- 
triguing triple point in the phase diagram of physical theories. Indeed, stick- 
ing to the existing literature (and refraining from dramatic statements on 
the historical evolution of theoretical physics) it appears that the subject of 
large-scale magnetization thrives and prosper at the crossroad of astrophysics, 
cosmology and theoretical high-energy physics. 

Following the kind invitation of Jnan Maharana and Maurizio Gasperini, I 
am delighted to contribute to this set of lectures whose guideline is dictated by 
the inspiring efforts of Gabriele Veneziano in understanding the fundamental 
forces of Nature. My voice joins the choir of gratitude proceeding from the 
whole physics community for the novel and intriguing results obtained by 
Gabriele through the various stages of his manifold activity. I finally ought 
to convey my personal thankfulness for the teachings, advices and generous 
clues received during the last fifteen years. 
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1.2 Length scales 

The typical magnetic field strengths, in the Universe, range from few fiG (in 
the case of galaxies and clusters), to few G (in the case of planets, like the 
earth or Jupiter) and up to 10 12 G in neutron stars. Magnetic fields are not 
only observed in planets and stars but also in the interstellar medium, in the 
intergalactic medium and, last but not least, in the intra-cluster medium. 

Magnetic fields whose correlation length is larger than the astronomical 
unit ( 1 AU = 1.49 x 10 13 cm) will be named large-scale magnetic fields. In fact, 
magnetic fields with approximate correlation scale comparable with the earth- 
sun distance are not observed (on the contrary, both the magnetic field of the 
sun and the one of the earth have a clearly distinguishable localized structure) . 
Moreover, in magnetohydrodynamics (MHD), the magnetic diffusivity scale 
(i.e. the scale below which magnetic fields are diffused because of the finite 
value of the conductivity) turns out to be, amusingly enough, of the order of 
the AU. 

1.3 The early history 

In the forties large-scale magnetic field had no empirical evidence. For in- 
stance, there was no evidence of magnetic fields associated with the galaxy 
as a whole with a rough correlation scale of 30kpc. More specifically, the 
theoretical situation can be summarized as follows. The seminal contributions 
of H. Alfven pQ convinced the community that magnetic fields can have a 
very large life-time in a highly conducting plasma. Later on, in the seventies, 
Alfven will be awarded by the Nobel prize "for fundamental work and discov- 
eries in magnetohydrodynamics with fruitful applications in different parts of 
plasma physics" . 

Using the discoveries of Alfven, Fermi [2] postulated, in 1949, the exis- 
tence of a large-scale magnetic field permeating the galaxy with approximate 
intensity of n G and, hence, in equilibrium with the cosmic rays[j 

Alfven [3] did not react positively to the proposal of Fermi, insisting, in a 
somehow opposite perspective, that cosmic rays are in equilibrium with stars 
and disregarding completely the possibility of a galactic magnetic field. Today 
we do know that this may be the case for low-energy cosmic rays but certainly 



1 Recall that lkpc = 3.085 x 10 21 cm. Moreover, IMpc = 10 3 kpc. The present size 
of the Hubble radius is H^ 1 = 1.2 x 10 28 cm = 4.1 x 10 3 Mpc for h = 0.73. 

2 In this contribution magnetic fields will be expressed in Gauss. In the SI units 
IT = 10 4 G. For practical reasons, in cosmic ray physics and in cosmology it is also 
useful to express the magnetic field in GeV 2 (in units h = c = 1). Recalling that 
the Bohr magneton is about 5.7 x 10~ 11 MeV/T the conversion factor will then be 
1G = 1.95 x 10" 20 GeV 2 . The use of Gauss (G) instead of Tesla (T) is justified by 
the existing astrophysical literature where magnetic fields are typically expressed 
in Gauss. 
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not for the most energetic ones around, and beyond, the knee in the cosmic 
ray spectrum. 

At the historical level it is amusing to notice that the mentioned contro- 
versy can be fully understood from the issue 75 of Physical Review where it is 
possible to consult the paper of Fermi [2] , the paper of Alfven [3] and even a 
paper by R. D. Richtmyer and E. Teller [4] supporting the views and doubts 
of Alfven. 

In 1949 Hiltner [5J and, independently, Hall [6] observed polarization of 
starlight which was later on interpreted by Davis and Greenstein [7] as an 
effect of galactic magnetic field aligning the dust grains. 

According to the presented chain of events it is legitimate to conclude that 

• the discoveries of Alfven were essential in the Fermi proposal who was pon- 
dering on the origin of cosmic rays in 1938 before leaving Italy because 
of the infamous fascist legislation; 

• the idea that cosmic rays are in equilibrium with the galactic magnetic 
fields (and hence that the galaxy possess a magnetic field) was essential in 
the correct interpretation of the first, fragile, optical evidence of galactic 
magnetization. 

The origin of the galactic magnetization, according to [5] , had to be somehow 
primordial. It should be noticed, for sake of completeness, that the observa- 
tions of Hiltner [5] and Hall [6] took place from November 1948 to January 
1949. The paper of Fermi [2] was submitted in January 1949 but it contains 
no reference to the work of Hiltner and Hall. This indicates the Fermi was 
probably not aware of these optical measurements. 

The idea that large-scale magnetization should somehow be the remnant 
of the initial conditions of the gravitational collapse of the protogalaxy idea 
was further pursued by Fermi in collaboration with S. Chandrasekar [8j [9] 
who tried, rather ambitiously, to connect the magnetic field of the galaxy to 
its angular momentum. 

1.4 The middle ages 

In the fifties various observations on polarization of Crab nebula suggested 
that the Milky Way is not the only magnetized structure in the sky. The 
effective new twist in the observations of large-scale magnetic fields was 
the development (through the fifties and sixties) of radio-astronomical tech- 
niques. From these measurements, the first unambiguous evidence of radio- 
polarization from the Milky Way (MW) was obtained (see [T0| and references 
therein for an account of these developments). 

It was also soon realized that the radio-Zeeman effect (counterpart of the 
optical Zeeman splitting employed to determine the magnetic field of the sun) 



3 The author is indebted with Prof. G. Cocconi who was so kind to share his 
personal recollections of the scientific discussions with E. Fermi. 
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could offer accurate determination of (locally very strong) magnetic fields 
in the galaxy. The observation of Lyne and Smith [TT] that pulsars could 
be used to determine the column density of electrons along the line of sight 
opened the possibility of using not only synchrotron emission as a diagnostic 
of the presence of a large-scale magnetic field, but also Faraday rotation. For 
a masterly written introduction to pulsar physics the reader may consult the 
book of Lyne and Smith [12] . 

In the seventies all the basic experimental tools for the analysis of galactic 
and extra-galactic magnetic fields were ready. Around this epoch also exten- 
sive reviews on the experimental endeavors started appearing and a very nice 
account could be found, for instance, in the review of Heiles [13]. 

It became gradually evident in the early eighties, that measurements of 
large-scale magnetic fields in the MW and in the external galaxies are two 
complementary aspects of the same problem. While MW studies can provide 
valuable informations concerning the local structure of the galactic magnetic 
field, the observation of external galaxies provides the only viable tool for the 
reconstruction of the global features of the galactic magnetic fields. 

Since the early seventies, some relevant attention has been paid not only 
to the magnetic fields of the galaxies but also to the magnetic fields of the 
clusters. A cluster is a gravitationally bound system of galaxies. The local 
group (i.e. our cluster containing the MW, Andromeda together with other 
fifty galaxies) is an irregular cluster in the sense that it contains fewer galaxies 
than typical clusters in the Universe. Other clusters (like Coma, Virgo) are 
more typical and are then called regular or Abell clusters. As an order of 
magnitude estimate, Abell clusters can contain 10 3 galaxies. 

1.5 New twists 

In the nineties magnetic fields have been measured in single Abell clusters but 
around the turn of the century these estimates became more reliable thanks 
to improved experimental techniques. In order to estimate magnetic fields in 
clusters, an independent knowledge of the electron density along the line of 
sight is needed. Recently Faraday rotation measurements obtained by radio 
telescopes (like VLA @) have been combined with independent measurements 
of the electron density in the intra-cluster medium. This was made possible by 
the maps of the x-ray sky obtained with satellites measurements (in particular 
ROSAT0). This improvement in the experimental capabilities seems to have 
partially settled the issue confirming the measurements of the early nineties 
and implying that also clusters are endowed with a magnetic field of [j,G 
strength which is not associated with individual galaxies [151 116] . 

4 The Very Large Array Telescope, consists of 27 parabolic antennas spread over a 
surface of 20 km 2 in Socorro (New Mexico) 

5 The ROegten SATellite (flying from June 1991 to February 1999) provided maps 
of the x-ray sky in the range 0.1-2.5 keV. A catalog of x-ray bright Abell clusters 
was compiled. 
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While entering the new millennium the capabilities of the observers are 
really confronted with a new challenge: the possibility that also superclus- 
ters are endowed with their own magnetic field. Superclusters are (loosely) 
gravitationally bound systems of clusters. An example is the local superclus- 
ter formed by the local group and by the VIRGO cluster. Recently a large 
new sample of Faraday rotation measures of polarized extragalactic sources 
has been compared with galaxy counts in Hercules and Perseus-Pisces (two 
nearby superclusters) [17] . First attempts to detect magnetic fileds associated 
with superclusters have been reported [TB]. A cautious and conservative ap- 
proach suggests that these fragile evidences must be corroborated with more 
conclusive observations (especially in the light of the, sometimes dubious, in- 
dependent determination of the electron density 0). However it is not excluded 
that as the nineties gave us a firmer evidence of cluster magnetism, the new 
millennium may give us more solid understanding of supercluster magnetism. 
In the present historical introduction various experimental techniques have 
been swiftly mentioned. A more extensive introductory description of these 
techniques can be found in 19J. 

1.6 Hopes for the future 

The hope for the near future is connected with the possibility of a next gener- 
ation radio-telescope. Along this line the SKA (Square Kilometer Array) has 
been proposed [16] (see also [20]). While the technical features of the instru- 
ment cannot be thoroughly discussed in the present contribution, it suffices 
to notice that the collecting area of the instrument, as the name suggest, will 
be of 10 6 m 2 . The specifications for the SKA require an angular resolution of 
0.1 arcsec at 1.4 GHz, a frequency capability of 0.1-25 GHz, and a field of 
view of at least 1 deg 2 at 1.4 GHz [5D]. The number of independent beams is 
expected to be larger than 4 and the number of instantaneous pencil beams 
will be roughly 100 with a maximum primary beam separation of about 100 
deg at low frequencies (becoming 1 deg at high frequencies, i.e. of the order of 
1 GHz). These specifications will probably allow full sky surveys of Faraday 
Rotation. 

The frequency range of SKA is rather suggestive if we compare it with the 
one of the Planck experiment [21j . Planck will operate in 9 frequency channels 
from 30 to, approximately, 900 GHz. While the three low-frequency channels 
(from 30 to 70 GHz) are not sensitive to polarization the six high-frequency 
channels (between 100 and 857 GHZ) will be definitely sensitive to CMB po- 
larization. Now, it should be appreciated that the Faraday rotation signal 



In [14] it was cleverly argued that informations on the plasma densities from 
direct observations can be gleaned from detailed multifrequency observations of 
few giant radio-galaxies (GRG) having dimensions up to 4 Mpc. The estimates 
based on this observation suggest column densities of electrons between 1CP 6 and 
1(T 5 cm" 3 . 
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decreases with the frequency v as v~ 2 . Therefore, for lower frequencies the 
Faraday Rotation signal will be larger than in the six high-frequency chan- 
nels. Consequently it is legitimate to hope for a fruitful interplay between the 
next generation of SKA-like radio-telescopes and CMB satellites. Indeed, as 
suggested above, the upper branch of the frequency capability of SKA almost 
overlaps with the lower frequency of Planck so that possible effects of large- 
scale magnetic fields on CMB polarization could be, with some luck, addressed 
with the combined action of both instruments. In fact, the same mechanism 
leading to the Faraday rotation in the radio leads to a Faraday rotation of the 
CMB provided the CMB is linearly polarized. These considerations suggest, as 
emphasized in a recent topical review, that CMB anisotropies are germane to 
several aspects of large-scale magnetization [22] . The considerations reported 
so far suggest that during the next decade the destiny of radio-astronomy and 
CMB physics will probably be linked together and not only for reasons of 
convenience. 

1.7 Few burning questions 

In this general and panoramic view of the history of the subject we started 
from the relatively old controversy opposing E. Fermi to H. Alfven with the 
still uncertain but foreseeable future developments. While the nature of the 
future developments is inextricably connected with the advent of new instru- 
mental capabilities, it is legitimate to remark that, in more than fifty years, 
magnetic fields have been detected over scales that are progressively larger. 
From the historical development of the subject a series of questions arises 
naturally: 

• what is the origin of large-scale magnetic fields? 

• are magnetic fields primordial as assumed by Fermi more than fifty years 
ago? 

• even assuming that large-scale magnetic fields are primordial, is there a 
theory for their generation? 

• is there a way to understand if large-scale magnetic fields are really pri- 
mordial? 

In what follows we will not give definite answers to these important questions 
but we shall be content of outlining possible avenues of new developments. 

The plan of the present lecture will be the following. In Sect. [2] the main 
theoretical problems connected with the origin of large-scale magnetic fields 
will be discussed. In Sect. [3] the attention will be focused on the problem of 
large-scale magnetic field generation in the framework of string cosmological 
model, a subject where the pre-big bang model, in its various incarnations, 
plays a crucial role. But, finally, large-scale magnetic fields are really primor- 
dial? Were they really present prior to matter-radiation equality? A modest 
approach to these important questions suggests to study the physics of mag- 
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netized CMB anisotropies which will be introduced, in its essential lines, in 
Sect. |U The concluding remarks are collected in Sect. 03 

2 Magnetogenesis 

While in the previous Section the approach has been purely historical, the 
experimental analysis of large-scale magnetic fields prompts a collection of 
interesting theoretical problems. They can be summarized by the following 
chain of evidences (see also [T5]): 

• In spiral galaxies magnetic fields follow the orientation of the spiral arms, 
where matter is clustered because of differential rotation. While there may 
be an asymmetry in the intensities of the magnetic field in the northern 
and southern emisphere (like it happens in the case of the Milky Way) the 
typical strength is in the range of the /i G. 

• Locally magnetic fields may even be in the mG range and, in this case, 
they may be detected through Zeeman splitting techniques. 

• In spiral galaxies the magnetic field is predominantly toroidal with a 
poloidal component present around the nucleus of the galaxy and extend- 
ing for, roughly, 100 pc. 

• The correlation scale of the magnetic field in spirals is of the order of 30 
kpc. 

• In elliptical galaxies magnetic fields have been measured at the fj, G level 
but the correlation scale is shorter than in the case of spirals: this is due 
to the different evolutionary history of elliptical galaxies and to their lack 
of differential rotation; 

• Abell clusters of galaxies exhibit magnetic fields present in the so-called 
intra-cluster medium: these fields, always at the fi G level, are not associ- 
ated with individual galaxies; 

• superclusters might also be magnetized even if, at the moment, conclusions 
are premature, as partially explained in Section Q] (see also [T5] and [H)]). 

The statements collected above rest on various detection techniques rang- 
ing from Faraday rotation, to synchrotron emission, to Zeeman splitting of 
clouds of molecules with an unpaired electron spin. The experimental evi- 
dence swiftly summarized above seems to suggest that different and distant 
objects have magnetic fields of comparable strength. The second suggestion 
seems also to be that the strength of the magnetic fields is, in the first (sim- 
plistic) approximation, independent on the physical scale. 

These empirical coincidences reminds a bit of one of the motivations of the 
standard hot big-bang model, namely the observation that the light elements 
are equally abundant in rather different parts of our Universe. The approxi- 
mate equality of the abundances implies that, unlike the heavier elements, the 
light elements have primordial origin. The four light isotopes D, 3 He, 4 He and 
7 Li are mainly produced at a specific stage of the hot big bang model named 
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nucleosynthesis occurring below the a typical temperature of 0.8 MeV when 
neutrinos decouple from the plasma and the neutron abundance evolves via 
free neutron decay [53]. The abundances calculated in the simplest big-bang 
nucleosythesis model agree fairly well with the astronomical observations. 

In similar terms it is plausible to argue that large-scale magnetic fields 
have comparable strengths at large scales because the initial conditions for 
their evolutions were the same, for instance at the time of the gravitational 
collapse of the protogalaxy. The way the initial conditions for the evolution of 
large-scale magnetic fields are set is generically named magnetogenesis |19j . 

There is another comparison which might be useful. Back in the seventies 
the so-called Harrison- Zeldovich spectrum was postulated. Later, with the de- 
velopments of inflationary cosmology the origin of a flat spectrum of curvature 
and density profiles has been justified on the basis of a period of quasi-de Sitter 
expansion named inflation. It is plausible that in some inflationary models not 
only the fluctuations of the geometry are amplified but also the fluctuations of 
the gauge fields. This happens if, for instance, gauge couplings are effectively 
dynamical. As the Harrison-Zeldovich spectrum can be used as initial condi- 
tion for the subsequent Newtonian evolution, the primordial spectrum of the 
gauge fields can be used as initial condition for the subsequent MHD evolu- 
tion which may lead, eventually, to the observed large-scale magnetic fields. 
The plan of the present section is the following. In Subsect. 12. II some general 
ideas of plasma physics will be summarized with particular attention to those 
tools that will be more relevant for the purposes of this lecture. In Subsect. 
12.21 the concept of dynamo amplification will be introduced in a simplified 
perspective. In Subsect. 12.31 it will be argued that the dynamo amplification, 
in one of its potential incarnations, necessitates some initial conditions or as 
we say in the jargon, some seed field. In Subsect. 12.41 a panoramic view of 
astrophysical seeds will be presented with the aim of stressing the common 
aspects of, sometimes diverse, physical mechanisms. Subsect. 12.51 and 12.61 the 
two basic approaches to cosmological magnetogenesis will be illustrated. In 
the first case (see Subsect. 12. 5|) magnetic fields are produced inside the Hub- 
ble radius at a given stage in the life of the Universe. In the second case (see 
Subsect. 12. 6|) vacuum fluctuations of the hypercharge field are amplified dur- 
ing an inflationary stage of expansion. Subsection 12.71 deals with the major 
problem of inflationary magnetogenesis, namely conformal (Weyl) invariance 
whose breaking will be one of the themes of string cosmological mechanisms 
for the generation of large-scale magnetic fields. 

2.1 Magnetized plasmas 

Large-scale magnetic fields evolve in a plasma, i.e. a system often illustrated 
as the fourth state of matter. As we can walk in the phase diagram of a given 
chemical element by going from the solid to the liquid and to the gaseous 
state with a series of diverse phase transitions, a plasma can be obtained by 
ionizing a gas. A typical example of weakly coupled plasma is therefore an 



Magnetic fields, strings and cosmology 



9 



ionized gas. Examples of strongly coupled plasmas can be found also in solid 
state physics. An essential physical scale that has to be introduced in the 
description of plasma properties is the so-called Debye length that will be 
discussed in the following paragraph. 

Different descriptions of a plasma exist and they range from effective fluid 
models of charged particles [211 E3 [26l [27] to kinetic approaches like the ones 
pioneered by Vlasov [28] and Landau [29]. From a physical point of view, a 
plasma is a system of charged particles which is globally neutral for typical 
length-scales larger than the Debye length Ad: 



where To is the kinetic temperature and no the mean charge density of the 
electron- ion system, i.e. n c ~ rij = no. For a test particle the Coulomb poten- 
tial will then have the usual Coulomb form but it will be suppressed, at large 
distances by a Yukawa term, i.e. e~ r / x ° . In the interstellar medium there are 
three kinds of regions which are conventionally defined: 

• H2 regions, where the Hydrogen is predominantly in molecular form (also 
denoted by HII); 

• H° regions (where Hydrogen is in atomic form); 

• and H + regions, where Hydrogen is ionized, (also denoted by HI). 

In the H + regions the typical temperature To is of the order of 10-20 eV while 
for no let us take, for instance, no ~ 3 x 10~ 2 cm~ 3 . Then Ad ~ 30km. 

For r ^> Ad the Coulomb potential is screened by the global effect of the 
other particles in the plasma. Suppose now that particles exchange momentum 
through two-body interactions. Their cross section will be of the order of 
O!g m /To and the mean free path will be £ m f p ~ To/(ag m no), i.e. recalling Eq. 
JI) Ad < ^mfp- This means that the plasma is a weakly collisional system 
which is, in general, not in local thermodynamical equilibrium and this is 
the reason why we introduced To as the kinetic (rather than thermodynamic) 
temperature. 

The last observation can be made even more explicit by defining another 
important scale, namely the plasma frequency which, in the system under 
discussion, is given by 



where m G is the electron mass. Notice that, in the interstellar medium (i.e. 
for no ~ 10~ 2 cm~ 3 ) Eq. gives a plasma frequency in the GHz range. This 
observation is important, for instance, in the treatment of Faraday rotation 
since the plasma frequency is typically much larger than the Larmor frequency 
i.e. 




(1) 
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"Be=^ 18.08 (^) kHz, (3) 

implying, for Bg ~ /xG, wb c — 20Hz. The same hierarchy holds also when the 
(free) electron density is much larger than in the interstellar medium, and, for 
instance, at the last scattering between electrons and photons for a rcdshift 
z dcc ~ 1100 (see Sect.0}. 

The plasma frequency is the oscillation frequency of the electrons when 
they are displaced from their equilibrium configuration in a background of ap- 
proximately fixed ions. Recalling that Vthcr — ^/To/m e is the thermal velocity 
of the charge carriers, the collision frequency uj c ~ fthcr/^mfp is always much 
smaller than u> pc ~ wthor/^D- Thus, in the idealized system described so far, 
the following hierarchy of scales holds: 

Ad < 4if P , w c < cjpc, (4) 

which means that before doing one collision the system undergoes many oscil- 
lations, or, in other words, that the mean free path is not the shortest scale in 
the problem. Usually one defines also the plasma parameter Af — »Iq , i.e. 
the number of particles in the Dcbye sphere. In the approximation of weakly 
coupled plasma Af -C 1 which also imply that the mean kinetic energy of the 
particles is larger than the mean inter-particle potential. 

The spectrum of plasma excitations is a rather vast subject and it will not 
strictly necessary for the following considerations (for further details see [24], 
[551 [55]). It is sufficient to remark that we can envisage, broadly speaking, two 
regimes that are physically different: 

• typical length-scales much larger than Ad and typical frequencies much 
smaller than io pc ; 

• typical length-scales smaller (or comparable) with Ad and typical frequen- 
cies much larger than uj pc . 

In the first situation reported above it can be shown that a single fluid de- 
scription suffices. The single fluid description is justified, in particular, for the 
analysis of the dynamo instability which occurs for dynamical times of the 
order of the age of the galaxy and length-scales larger than the kpc. In the 
opposite regime, i.e. u > u> pc and L > Ad the single fluid approach breaks 
down and a multi-fluid description is mandatory This is, for instance, the 
branch of the spectrum of plasma excitation where the displacement current 
(and the related electromagnetic propagation) cannot be neglected. A more 
reliable description is provided, in this regime, by the Vlasov-Landau (i.e. 
kinetic) approach [551 US] (see also [55]). 

Consider, therefore, a two-fluid system of electrons and protons. This sys- 
tem will be described by the continuity equations of the density of particles, 
i.e. 

^ + V • (n e v e ) = 0, ^ + V • (n p v p ) = 0, (5) 
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and by the momentum conservation equations 





' d 






m c n c 


dt ' 




v c = -en c 



E + v c x B 



E + v p x B 



- Vp c - C cp , 

- Vp p - c pc . 



(6) 
(7) 



Equations ([5| , © and (0 must be supplemented by Maxwell equations read- 
ing, in this case 



V • E = 47re(n p — n c ), 

V B = 0, 

dB 

VxE+ — =0, 
ot 

dE 

V x B = — + 4ne(n p v p - n c v c ). 



(8) 
(9) 

(10) 
(11) 



The two fluid system of equations is rather useful to discuss various phe- 
nomena like the propagation of electromagnetic excitations at finite charge 
density both in the presence and in the absence of a background magnetic 
field [2H [26]. The previous observation implies that a two-fluid treat- 
ment is mandatory for the description of Faraday rotation of the Cosmic Mi- 
crowave Background (CMB) polarization. This subject will not be specifically 
discussed in the present lecture (see, for further details, [30] and references 
therein) . 

Instead of treating the two fluids as separated, the plasma may be consid- 
ered as a single fluid defined by an appropriate set of global variables: 



J = e(n p v p - n c v c ), 
P q = e(n p - n e ), 
Pro = (m c n c + m p n p ), 
m c n c v c + n p m p v p 
m c n c + m p n p 



v = 



(12) 
(13) 
(14) 

(15) 



where J is the global current and p q is the global charge density; /? m is the 
total mass density and v is the so-called bulk velocity of the plasma. From the 
definition of the bulk velocity it is clear that v is the centre-of-mass velocity 
of the electron-ion system. The interesting case is the one where the plasma 
is globally neutral, i.e. n c ~ n p — no, implying, from Maxwell and continuity 
equations the following equations 



V-£ = 0, 



V • J = 0, 



V • B = 0. 



(16) 



The equations reported in Eq. (fl^]) are the first characterization of MHD 
equations, i.e. a system where the total current as well as the electric and 
magnetic fields are all solenoidal. The remaining equations allow to obtain 
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the relevant set of conditions describing the long wavelength modes of the 
magnetic field i.e. 

VxB = 4ttJ, (17) 

vx *; = --. (is) 

In Eq. ([T7| , the contribution of the displacement current has been neglected 
for consistency with the solenoidal nature of the total current (see Eq. (TlGl) ). 
Two other relevant equations can be obtained by summing and subtracting 
the momentum conservation equations, i.e. Eqs. J6]) and ([7]). The result of this 
procedure is 



dv 

— + v ■ Vv 

at 



= J x B - VP (19) 



E + v x B = - + — (J x B - Vp e ), (20) 
(j en q 

where n q ~ ti ~ n c and P — p c + p p . Equation (| 19[) is derived from the 
sum of Eqs. |6]) and ([7|) and in (fT9|) J x £? is the Lorentz force term which is 
quadratic in the magnetic field. In fact using Eq. (|17j) 

J x B — —(V x B) x B. (21) 

Note that to derive Eq. (|2l7|) the limit m e /m p must be taken, at some 
point. There are some caveats related to this procedure since viscous and 
collisional effects may be relevant [25] . Equation (|20|) is sometimes called one- 
fluid generalized Ohm law. In Eq. (|2H)l the term J x B is nothing but the Hall 
current and Vp e is often called thermoelectric term. Finally the term J /a is 
the resistivity term and a is the conductivity of the one-fluid description. In 
Eq. (|20[) the pressure has been taken to be isotropic. Neglecting, the Hall and 
thermoelectric terms (that may play, however, a role in the Biermann battery 
mechanism for magnetic field generation) the Ohm law takes the form 

J = (j(E + v x B). (22) 



Using Eq. (|22|) together with Eq. (| 1 T[) it is easy to show that the Ohmic electric 
field is given by 

E = V * B -v x B. (23) 

47TCT 

Using then Eq. (I23p into Eq. (fT5|) and exploiting known vector identities we 
can get the canonical form of the magnetic diffusivity equation 

OB 1 

- = V x (* x B) + — V 2 B, (24) 

which is the equation to be used to discuss the general features of the dynamo 
instability. 

MHD can be studied into two different (but complementary) limits 
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• the ideal (or superconducting) limit where the conductivity is set to infinity 
(i.e. the a — > oo limit); 

• the real (or resistive) limit where the conductivity is finite. 

The plasma description following from MHD can be also phrased in terms 
of the conservation of two interesting quantities, i.e. the magnetic flux and 
the magnetic helicity [27J, [31] : 

In Eq. (|25p . S is an arbitrary closed surface that moves with the plasma. In 
the ideal MHD limit the magnetic flux is exactly conserved and the the flux 
is sometimes said to be frozen into the plasma element. In the same limit also 
the magnetic helicity is conserved. In the resistive limit the magnetic flux and 
helicity arc dissipated with a rate proportional to 1/a which is small provided 
the conductivity is sufficiently high. The term appearing at the right hand 
side off Eq. (|2"6"|) is called magnetic gyrotropy. 

The conservation of the magnetic helicity is a statement on the conserva- 
tion of the topological properties of the magnetic flux lines. If the magnetic field 
is completely stochastic, the magnetic flux lines will be closed loops evolving 
independently in the plasma and the helicity will vanish. There could be, how- 
ever, more complicated topological situations where a single magnetic loop is 
twisted (like some kind of Mobius stripe) or the case where the magnetic loops 
are connected like the rings of a chain. In both cases the magnetic helicity 
will not be zero since it measures, essentially, the number of links and twists 
in the magnetic flux lines. The conservation of the magnetic flux and of the 
magnetic helicity is a consequence of the fact that, in ideal MHD, the Ohmic 
electric field is always orthogonal both to the bulk velocity field and to the 
magnetic field. In the resistive MHD approximation this is no longer true [27] . 

2.2 Dynamos 

The dynamo theory has been developed starting from the early fifties through 
the eighties and various extensive presentations exist in the literature [32[ I33|, 
134) . Generally speaking a dynamo is a process where the kinetic energy of the 
plasma is transferred to magnetic energy. There arc different sorts of dynamos. 
Some of the dynamos that are currently addressed in the existing literature are 
large-scale dynamos, small-scale dynamos, nonlinear dynamos, a-dynamos... 

It would be difficult, in the present lecture, even to review such a vast lit- 
erature and, therefore, it is more appropriate to refer to some review articles 
where the modern developments in dynamo theory and in mean field elec- 
trodynamics are reported [35[ [37] . As a qualitative example of the dynamo 



1 

47T<7 
1 
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V x V x B ■ dU, 
d 3 xB ■ V x B. 



(25) 
(26) 
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action it is practical do discuss the magnetic diffusivity equation obtained, 
from general considerations, in Eq. (I24|) . 

Equation (pM)) simply stipulates that the first time derivative of the mag- 
netic fields intensity results from the balance of two (physically different) 
contributions. The first term at the right hand side of Eq. ([24]) is the the 
dynamo term and it contains the bulk velocity of the plasma v. If this term 
dominates the magnetic field may be amplified thanks to the differential ro- 
tation of the plasma. The dynamo term provides then the coupling allowing 
the transfer of the kinetic energy into magnetic energy. The second term at 
the right hand side of Eq. (|24[) is the magnetic diffusivity whose effect is to 
damp the magnetic field intensity. Defining then as L the typical scale of spa- 
tial variation of the magnetic field intensity, the typical time scale of resistive 
phenomena turns out to be 

t a ~ 4naL 2 . (27) 

In a non-relativistic plasma the conductivity a goes typically as T 3 / 2 [MISS]. 
In the case of planets, like the earth, one can wonder why a sizable magnetic 
field can still be present. One of the theories is that the dynamo term regen- 
erates continuously the magnetic field which is dissipated by the diffusivity 
term |32j . In the case of the galactic disk the value of the conductivity is 
given by a ~ 7 x 10 _7 Hz. Thus, for L ~ kpc t a ~ 10 9 (L/kpc) 2 sec. 

Equation (|27|) can also give the typical resistive length scale once the time- 
scale of the system is specified. Suppose that the time-scale of the system is 
given by tjj ~ H^ 1 ~ 10 18 sec where Hq is the present order of magnitude of 
the Hubble parameter. Then 




leading to L a ~ AU. The scale ([2"5]) gives then the upper limit on the diffusion 
scale for a magnetic field whose lifetime is comparable with the age of the 
Universe at the present epoch. Magnetic fields with typical correlation scale 
larger than L a are not affected by resistivity. On the other hand, magnetic 
fields with typical correlation scale L < L a are diffused. The value L a ~ AU 
is consistent with the phenomenological evidence that there are no magnetic 
fields coherent over scales smaller than 10~ 5 pc. 

The dynamo term may be responsible for the origin of the magnetic field 
of the galaxy. The galaxy has a typical rotation period of 3 x 10 8 yrs and 
comparing this figure with the typical age of the galaxy, C(10 10 yrs), it can 
be appreciated that the galaxy performed about 30 rotations since the time 
of the protogalactic collapse. 

The effectiveness of the dynamo action depends on the physical properties 
of the bulk velocity field. In particular, a necessary requirement to have a 



7 It is common use in the astrophysical applications to work directly with r\ — 
(Ana) -1 . In the case of the galactic disks n = 10 26 cm 2 Hz. 
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potentially successful dynamo action is that the velocity field is non-mirror- 
symmetric or that, in other words, (v ■ V x v) ^ 0. Let us see how this 
statement can be made reasonable in the framework of Eq. (pM)) . From Eq. 
(|24[) the usual structure of the dynamo term may be derived by carefully 
averaging over the velocity filed according to the procedure of [4TJ [42] . By 
assuming that the motion of the fluid is random and with zero mean velocity 
the average is taken over the ensemble of the possible velocity fields. In more 
physical terms this averaging procedure of Eq. (f24"|) is equivalent to average 
over scales and times exceeding the characteristic correlation scale and time 
To of the velocity field. This procedure assumes that the correlation scale of 
the magnetic field is much bigger than the correlation scale of the velocity 
field which is required to be divergence- less (V ■ v = 0). In this approximation 
the magnetic diffusivity equation can be written as: 

f) F? 1 

-=a(VxB) + —^B, (29) 

where 

a = ~(v ■ V x v), (30) 

is the so-called a-term in the absence of vorticity. In Eqs. (f29)) - (|3"0")) B is 
the magnetic field averaged over times longer that to which is the typical 
correlation time of the velocity field. 

The fact that the velocity field must be globally non-mirror symmetric 
[3"5] suggests, already at this qualitative level, the deep connection between 
dynamo action and fully developed turbulence. In fact, if the system would 
be, globally, invariant under parity transformations, then, the a term would 
simply be vanishing. This observation may also be related to the turbulent 
features of cosmic systems. In cosmic turbulence the systems are usually ro- 
tating and, moreover, they possess a gradient in the matter density (think, for 
instance, to the case of the galaxy). It is then plausible that parity is broken 
at the level of the galaxy since terms like V/? m • V x v are not vanishing 33J. 

The dynamo term, as it appears in Eq. (|29|) , has a simple electrodynamical 
meaning, namely, it can be interpreted as a mean ohmic current directed along 
the magnetic field : 

J = aB. (31) 

Equation stipulates that an ensemble of screw-like vortices with zero mean 
helicity is able to generate loops in the magnetic flux tubes in a plane orthog- 
onal to the one of the original field. As a simple (and known) application of 
Eq. (I29p , it is appropriate to consider the case where the magnetic field profile 
is given by a sort of Chern-Simons wave 

B x (z,t) =f(t)sixikz, B v = /(f)cosfcz, B z (k,t) = 0. (32) 

For this profile the magnetic gyrotropy is non-vanishing, i.e. B ■ V x B = 
kf 2 (t). From Eq. (|2"9"|) , using Eq. (p{2"|) f(t) obeys the following equation 
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dt \ Attct J 



f (33) 



admits exponentially growing solutions for sufficiently large scales, i.e. k < 
4n\a\a. Notice that in this naive example the a term is assumed to be con- 
stant. However, as the amplification proceeds, a may develop a dependence 
upon \B\ 2 , i.e. a -> a (l - £|-B| 2 )a [l - £/ 2 (*)]. In the case of Eq. ||3SJ) this 
modification will introduce non-linear terms whose effect will be to stop the 
growth of the magnetic field. This regime is often called saturation of the dy- 
namo and the non-linear equations appearing in this context are sometimes 
called Landau equations [33] in analogy with the Landau equations appearing 
in hydrodynamical turbulence. 

In spite of the fact that in the previous example the velocity field has 
been averaged, its evolution obeys the Navier-Stokes equation which we have 
already written but without the diffusion term 



— + (v ■ V)i> - vV z v 
at 



-VP + JxB, (34) 



where v is the thermal viscosity coefficient. There are idealized cases where the 
Lorentz force term can be neglected. This is the so-called force free approxi- 
mation. Defining the kinetic helicity as fl = V x v, the magnetic diffusivity 
and Navier-Stokes equations can be written in a rather simple and symmetric 
form 

OB 1 
-=Vx(vxB) + — V**, 

Ofl o 

— = V x (v x fl) + vV 2 fl. (35) 

In MHD various dimensionless ratios can be defined. The most frequently 
used are the magnetic Reynolds number, the kinetic Reynolds number and 
the Prandtl number: 

R m = vL B <j, (36) 
R=^, (37) 

*- (38) 

where Lb and L v are the typical scales of variation of the magnetic and veloc- 
ity fields. If R m 3> 1 the system is said to be magnetically turbulent. If R 3> 1 
the system is said to be kinetically turbulent. In realistic situations the plasma 
is both kinetically and magnetically turbulent and, therefore, the ratio of the 
two Reynolds numbers will tell which is the dominant source of turbulence. 
There have been, in recent years, various studies on the development of mag- 
netized turbulence (see, for instance, [27]) whose features differ slightly from 
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the ones of hydrodynamic turbulence. While the details of this discussion will 
be left aside, it is relevant to mention that, in the early Universe, turbulence 
may develop. In this situation a typical phenomenon, called inverse cascade, 
can take place. A direct cascade is a process where energy is transferred from 
large to small scales. Even more interesting, for the purposes of the present 
lecture, is the opposite process, namely the inverse cascade where the energy 
transfer goes from small to large length-scales. One can also generalize the 
the concept of energy cascade to the cascade of any conserved quantity in the 
plasma, like, for instance, the helicity. Thus, in general terms, the transfer 
process of a conserved quantity is a cascade. 

The concept of cascade (either direct or inverse) is related with the concept 
of turbulence, i.e. the class of phenomena taking place in fluids and plasmas at 
high Reynolds numbers. It is very difficult to reach, with terrestrial plasmas, 
the physical situation where the magnetic and the kinetic Reynolds numbers 
are both large but, in such a way that their ratio is also large i.e. 

R 

R m > 1, R > 1, Pr = > 1. (39) 

R 

The physical regime expressed through Eqs. (|39p rather common in the early 
Universe. Thus, MHD turbulence is probably one of the key aspects of mag- 
netized plasma dynamics at very high temperatures and densities. Consider, 
for instance, the plasma at the electroweak epoch when the temperature was 
of the order of 100 GeV. One can compute the Reynolds numbers and the 
Prandtl number from their definitions given in Eqs. (|36p - (|38p . In particular, 

R m -10 17 , R=10 n , Pr~ 10 6 , (40) 

which can be obtained from Eqs. ([36|) - (j38)) using as fiducial parameters v ~ 
0.1, crT/a, v ~ (aT)- 1 and L ~ 0.01 H~J ~ 0.03 cm for T ~ 100 GeV. 

If an inverse energy cascade takes place, many (energetic) magnetic do- 
mains coalesce giving rise to a magnetic domain of larger size but of smaller 
energy. This phenomenon can be viewed, in more quantitative terms, as an 
effective increase of the correlation scale of the magnetic field. This consider- 
ation plays a crucial role for the viability of mechanisms where the magnetic 
field is produced in the early Universe inside the Hubble radius (see Subsect. 

123). 

2.3 Initial conditions for dynamos 

According to the qualitative description of the dynamo instability presented 
in the previous subsection, the origin of large-scale magnetic fields in spiral 
galaxies can be reduced to the three keywords: seeding, amplification and 
ordering. The first stage, i.e. the seeding, is the most controversial one and 
will be briefly reviewed in the following sections of the present review. In more 
quantitative terms the amplification and the ordering may be summarized as 
follows: 
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• during the 30 rotations performed by the galaxy since the protogalactic 
collapse, the magnetic field should be amplified by about 30 e-folds; 

• if the large scale magnetic field of the galaxy is, today, O(^G) the magnetic 
field at the onset of galactic rotation might have been even 30 e-folds 
smaller, i.e. 0(1O _19 G) over a typical scale of 30-100 kpc; 

• assuming perfect flux freezing during the gravitational collapse of the pro- 
togalaxy (i.e. a — > oo) the magnetic field at the onset of gravitational 
collapse should be 0(1O~ 23 ) G over a typical scale of 1 Mpc. 

This picture is oversimplified and each of the three steps mentioned above can 
be questioned. In what follows the main sources of debate, emerged in the last 
ten years, will be briefly discussed. 

There is a simple way to relate the value of the magnetic fields right after 
gravitational collapse to the value of the magnetic field right before gravi- 
tational collapse. Since the gravitational collapse occurs at high conductivity 
the magnetic flux and the magnetic helicity are both conserved (see, in partic- 
ular, Eq. ([25]) ). Right before the formation of the galaxy a patch of matter of 
roughly 1 Mpc collapses by gravitational instability. Right before the collapse 
the mean energy density of the patch, stored in matter, is of the order of the 
critical density of the Universe. Right after collapse the mean matter density 
of the protogalaxy is, approximately, six orders of magnitude larger than the 
critical density. 

Since the physical size of the patch decreases from 1 Mpc to 30 kpc the 
magnetic field increases, because of flux conservation, of a factor (p a /pb) 2 ^ 3 ~ 
10 4 where p a and pb are, respectively the energy densities right after and right 
before gravitational collapse. The correct initial condition in order to turn on 
the dynamo instability would be \B\ ~ 10~ 23 Gauss over a scale of 1 Mpc, 
right before gravitational collapse. 

The estimates presented in the last paragraph are based on the (rather 
questionable) assumption that the amplification occurs over thirty e-folds 
while the magnetic flux is completely frozen in. In the real situation, the 
achievable amplification is much smaller. Typically a good seed would not be 
10 -19 G after collapse (as we assumed for the simplicity of the discussion) but 
rather [35] 




The galactic rotation period is of the order of 3 x 10 8 yrs. This scale 
should be compared with the typical age of the galaxy. All along this rather 
large dynamical time-scale the effort has been directed, from the fifties, to 
the justification that a substantial portion of the kinetic energy of the system 
(provided by the differential rotation) may be converted into magnetic energy 
amplifying, in this way, the seed field up to the observed value of the magnetic 
field, for instance in galaxies and in clusters. In recent years a lot of progress 
has been made both in the context of the small and large-scale dynamos [36J 
[37] (see also [38l[39l|40]). This progress was also driven by the higher resolution 
of the numerical simulations and by the improvement in the understanding of 




\B\ > 10~ 13 G. 



(41) 
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the largest magnetized system that is rather close to us, i.e. the sun [37]. More 
complete accounts of this progress can be found in the second paper of Ref. 
[36] and, more comprehensively, in Ref. [37] . Apart from the aspects involving 
solar physics and numerical analysis, better physical understanding of the role 
of the magnetic helicity in the dynamo action has been reached. This point 
is crucially connected with the two conservation laws arising in MHD, i.e. 
the magnetic flux and magnetic helicity conservations whose relevance has 
been already emphasized, respectively, in Eqs. ([25]) and (|26|) . Even if the rich 
interplay between small and large scale dynamos is rather important, let us 
focus on the problem of large-scale dynamo action that is, at least superficially, 
more central for the considerations developed in the present lecture. 

Already at a qualitative level it is clear that there is a clash between the 
absence of mirror-symmetry of the plasma, the quasi-exponential amplification 
of the seed and the conservation of magnetic flux and helicity in the high (or 
more precisely infinite) conductivity limit. The easiest clash to understand, 
intuitively, is the flux conservation versus the exponential amplification: both 
flux freezing and exponential amplification have to take place in the same 
superconductive (i.e. a^ 1 — > 0) limit. The clash between helicity conservation 
and dynamo action can be also understood in general terms: the dynamo 
action implies a topology change of the configuration since the magnetic flux 
lines cross each other constantly 36 . 

One of the recent progress in this framework is a more consistent formula- 
tion of the large-scale dynamo problem (36J, [37] : large scale dynamos produces 
small scale helical fields that quench (i.e. prematurely saturate) the a effect. 
In other words, the conservation of the magnetic helicity can be seen, accord- 
ing to the recent view, as a fundamental constraint on the dynamo action. 
In connection with the last point, it should be mentioned that, in the past, a 
rather different argument was suggested 46J: it was argued that the dynamo 
action not only leads to the amplification of the large-scale field but also of the 
random field component. The random field would then suppress strongly the 
dynamo action. According to the considerations based on the conservation of 
the magnetic helicity this argument seems to be incorrect since the increase 
of the random component would also entail and increase of the rate of the 
topology change, i.e. a magnetic helicity non-conservation. 

The possible applications of dynamo mechanism to clusters is still under 
debate and it seems more problematic. The typical scale of the gravitational 
collapse of a cluster is larger (roughly by one order of magnitude) than the 
scale of gravitational collapse of the protogalaxy. Furthermore, the mean mass 
density within the Abell radius ( ~ 1.5/i -1 Mpc) is roughly 10 3 larger than the 
critical density. Consequently, clusters rotate much less than galaxies. Recall 
that clusters are formed from peaks in the density field. The present overden- 
sity of clusters is of the order of 10 3 . Thus, in order to get the intra-cluster 
magnetic field, one could think that magnetic flux is exactly conserved and, 
then, from an intcrgalactic magnetic field \B\ > 10 -9 G an intra cluster mag- 
netic field |£?| > 10~ 7 G can be generated. This simple estimate shows why it 
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is rather important to improve the accuracy of magnetic field measurements 
in the intra-cluster medium: the change of a single order of magnitude in the 
estimated magnetic field may imply rather different conclusions for its origin. 



2.4 Astrophysical mechanisms 

Many (if not all) the astrophysical mechanisms proposed so far are related 
to what is called, in the jargon, a battery. In short, the idea is the following. 
The explicit form of the generalized Ohmic electric field in the presence of 
thermoelectric corrections can be written as in Eq. (|20ll where we set n q = n c 
to stick to the usual conventional 

E = -vxB + ^--^. (42) 
47rcr en e 



By comparing Eq. (|23|) with Eq. (1421) . it is clear that the additional term 
at the right hand side, receives contribution from a temperature gradient. 
In fact, restoring for a moment the Boltzmann constant ks we have that 
since P c = ^ n c T c . the additional term depends upon the gradients of the 
temperature, hence the name thermoelectric. It is interesting to see under 
which conditions the curl of the electric field receives contribution from the 
thermoelectric effect. Taking the curl of both sides of Eq. (j42|) we obtain 



V x E = + V(v x B) - X 2 VPc = -gg (43) 

Ana erig at 

where the second equality is a consequence of Maxwell's equations. From 
Eq. (|43|) it is clear that the evolution of the magnetic field inherits a source 
term iff the gradients in the pressure and electron density are not parallel. If 
VP C || Vn c a fully valid solution of Eq. (|43|) is B — 0. In the opposite case 
a seed magnetic field is naturally provided by the thermoelectric term. The 
usual (and rather general) observation that one can make in connection with 
the geometrical properties of the thermoelectric term is that cosmic ionization 
fronts may play an important role. For instance, when quasars emit ultraviolet 
photons, cosmic ionization fronts are produced. Then the intcrgalactic medium 
may be ionized. It should also be recalled, however, that the temperature 
gradients are usually normal to the ionization front. In spite of this, it is also 
plausible to think that density gradients can arise in arbitrary directions due 
to the stochastic nature of density fluctuations. 



For simplicity, we shall neglect the Hall contribution arising in the generalized 
Ohm law. The Hall contribution would produce, in Eq. (|42|l a term J x B/n c e 
that is of higher order in the magnetic field and that is proportional to the Lorentz 
force. The Hall term will play no role in the subsequent considerations. However, 
it should be borne in mind that the Hall contribution may be rather interesting 
in connection with the presence of strong magnetic fields like the ones of neutron 
stars (i.e. 10 13 G). This occurrence is even more interesting since in the outer 
regions of neutron stars strong density gradients are expected. 
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In one way or in another, astrophysical mechanisms for the generation of 
magnetic fields use an incarnation of the thermoelectric effect [13] (see also 
[HISS]). In the sixties and seventies, for instance, it was rather popular to 
think that the correct "geometrical" properties of the thermoelectric term may 
be provided by a large-scale vorticity. As it will also be discussed later, this 
assumption seems to be, at least naively, in contradiction with the formulation 
of inflationary models whose prediction would actually be that the large-scale 
vector modes are completely washed-out by the expansion of the Universe. 
Indeed, all along the eighties and nineties the idea of primordial vorticity 
received just a minor attention. 

The attention then focused on the possibility that objects of rather small 
size may provide intense seeds. After all we do know that these objects may 
exist. For instance the Crab nebula has a typical size of a roughly 1 pc and 
a magnetic field that is a fraction of the m G. These seeds will then combine 
and diffuse leading, ultimately, to a weaker seed but with large correlation 
scale. This aspect, may be, physically, a bit controversial since we do observe 
magnetic fields in galaxies and clusters that are ordered over very large length 
scales. It would then seem necessary that the seed fields produced in a small 
object (or in several small objects) undergo some type of dynamical self- 
organization whose final effect is a seed coherent over length-scales 4 or 5 
orders of magnitude larger than the correlation scale of the original battery. 

An interesting idea could be that qualitatively different batteries lead to 
some type of conspiracy that may produce a strong large scale seed. In [43] it 
has been suggested that Population III stars may become magnetized thanks 
to a battery operating at stellar scale. Then if these stars would explode as 
supernovae (or if they would eject a magnetized stellar wind) the pre-galactic 
environment may be magnetized and the remnants of the process incorporated 
in the galactic disc. In a complementary perspective, a similar chain of events 
may take place over a different physical scale. A battery could arise, in fact 
in active galactic nuclei at high red-shift. Then the magnetic field could be 
ejected leading to intense fields in the lobes of "young" radio-galaxies. These 
fields will be somehow inherited by the "older" disc galaxies and the final seed 
field may be, according to [43] as large as 10 -9 G at the pre-galactic stage. 

In summary we can therefore say that: 

• both the primordial and the astrophysical hypothesis for the origin of the 
seeds demand an efficient (large-scale) dynamo action; 

• due to the constraints arising from the conservation of magnetic helicity 
and magnetic flux the values of the required seed fields may turn out to be 
larger than previously thought at least in the case when the amplification 
is only driven by a large-scale dynamo action^; 



The situation may change if the magnetic fields originate from the combined 
action of small and large scale dynamos like in the case of the two-step process 
described in [43] . 
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• magnetic flux conservation during gravitational collapse of the protogalaxy 
may increase, by compressional amplification, the initial seed of even 4 
orders of magnitude; 

• compressional amplification, as well as large-scale dynamo, are much less 
effective in clusters: therefore, the magnetic field of clusters is probably 
connected to the specific way the dynamo saturates, and, in this sense, 
harder to predict from a specific value of the initial seed. 

2.5 Magnetogenesis: inside the Hubble radius 

One of the weaknesses of the astrophysical hypothesis is connected with the 
smallness of the correlation scale of the obtained magnetic fields. This type 
of impasse led the community to consider the option that the initial condi- 
tions for the MHD evolution are dictated not by astrophysics but rather by 
cosmology. The first ones to think about cosmology as a possible source of 
large-scale magnetization were Zeldovich [ITJ EE] ; an d Harrison [151 HOI H2] ■ 

The emphasis of these two authors was clearly different. While Zeldovich 
thought about a magnetic field which is uniform (i.e. homogeneous and ori- 
ented, for instance, along a specific Cartesian direction) Harrison somehow 
anticipated the more modern view by considering the possibility of an inho- 
mogeneous magnetic field. In the scenario of Zeldovich the uniform magnetic 
field would induce a slight anisotropy in the expansion rate along which the 
magnetic field is aligned. So, for instance, by considering a constant (and uni- 
form) magnetic field pointing along the x Cartesian axis, the induced geometry 
compatible with such a configuration will fall into the Bianchi-I class 

ds 2 = dt 2 - a 2 {t)dx 2 - b 2 {t)[dy 2 + dz 2 ]. (44) 

By solving Einstein equations in this background geometry it turns out that, 
during a radiation dominated epoch, the expansion rates along the x and the 
y — z plane change and their difference is proportional to the magnetic energy 
density [UJ EH] ■ This observation is not only relevant for magnetogenesis but 
also for Cosmic Microwave Background (CMB) anisotropies since the differ- 
ence in the expansion rate turns out to be proportional to the temperature 
anisotropy. While we will get back to this point later, in Section 21 as far as 
magnetization is concerned we can just remark that the idea of Zeldovich was 
that a uniform magnetic field would modify the initial condition of the stan- 
dard hot big bang model where the Universe would start its evolution already 
in a radiation-dominated phase. 

The model of Harrison [15J [SHI [H] is, in a sense, more dynamical. Fol- 
lowing earlier work of Bicrmann [52] . Harrison thought that inhomogencous 
MHD equations could be used to gennerate large-scale magnetic fields provided 
the velocity field was turbulent enough. The Biermann battery was simply a 
battery (as the ones described above in this session) but operating prior to 
decoupling of matter and radiation. The idea of Harrison was instead that 
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vorticity was already present so that the effective MHD equations will take 
the form 

^-(a 2 n + —B) = -^—^ 2 B, (45) 

OT m p 47T(TTOp 

where, as previously defined, 17 = V x v and m p is the ion mass. Equation 
([45)) is written in a conformally flat Friedmann-Robertson- Walker metric of 
the form 

ds 2 = G^dx^dx" = a 2 (r)[rfr 2 - dx\ (46) 

where r is the conformal time coordinate and where, in the conformally flat 
case, G^y = a 2 (r)?7 

y,v i "H^v being the four-dimensional Minkowski metric. If 
we now postulate that some vorticity was present prior to decoupling, then 
Eq. can be solved and the magnetic field can be related to the initial 
vorticity as 

B ~-!?W^) 2 . (47) 

If the estimate of the vorticity is made prior to equality (as originally 
suggested by Harrison [49]) of after decoupling as also suggested, a bit later, 
in Ref. [53], the result can change even by two orders of magnitude. Prior to 
equality \f2(i) ~ 0.1/f and, therefore, |-B oq | ~ 1CU 21 G. If a similar estimate is 
made after decoupling the typical value of the generated magnetic field is of 
the order of 10~ 18 G. So, in this context, the problem of the origin of magnetic 
fields is circumvented by postulating an appropriate form of vorticity whose 
origin must be explained. 

The Harrison mechanism is just one of the first examples of magnetic field 
generation inside the Hubble radius. In cosmology we define the Hubble radius 
as the inverse of the Hubble parameter, i.e. r# = H~ l (t). The first possibility 
we can think of implies that magnetic fields are produced, at a given epoch 
in the life of the Universe, inside the Hubble radius, for instance by a phase 
transition or by any other phenomenon able to generate a charge separation 
and, ultimately, an electric current. In this context, the correlation scale of 
the field is much smaller that the typical scale of the gravitational collapse of 
the proto-galaxy which is of the order of the Mpc. In fact, if the Universe is 
decelerating and if the correlation scale evolves as the scale factor, the Hubble 
radius grows much faster than the correlation scale. Of course, one might 
invoke the possibility that the correlation scale of the magnetic field evolves 
more rapidly than the scale factor. A well founded physical rationale for this 
occurrence is what is normally called inverse cascade, i.e. the possibility that 
magnetic (as well as kinetic) energy density is transferred from small to large 
scales. This implies, in real space, that (highly energetic) small scale magnetic 
domains may coalesce to form magnetic domains of smaller energy but over 
larger scales. In the best of all possible situations, i.e. when inverse cascade is 
very effective, it seems rather hard to justify a growth of the correlation scale 
that would eventually end up into a Mpc scale at the onset of gravitational 
collapse. 
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INSIDE THE HUBBLE RADIUS 




Fig. 1. Evolution of the correlation scale for magnetic fields produced inside the 
Hubble radius. The horizontal thick dashed line marks the end of the radiation- 
dominated phase and the onset of the matter-dominated phase. The horizontal thin 
dashed line marks the moment of e + -e~ annihilation (see also footnoote 2). The 
full (vertical) lines represent the evolution of the Hubble radius during the different 
stages of the life of the Universe. The dashed (vertical) lines illustrate the evolution 
of the correlation scale of the magnetic fields. In the absence of inverse cascade the 
evolution of the correlation scale is given by the (inner) vertical dashed lines. If 
inverse cascade takes place the evolution of the correlation scale is faster than the 
first power of the scale factor (for instance a 5 ^ 3 ) but always slower than the Hubble 
radius. 

In Fig. [1] we report a schematic illustration of the evolution of the Hubble 
radius Rh and of the correlation scale of the magnetic field as a function of 
the scale factor. In Fig. [1] the horizontal dashed line simply marks the end of 
the radiation-dominated phase and the onset of the matter dominated phase: 
while above the dashed line the Hubble radius evolves as a 2 (where a is the 
scale factor), below the dashed line the Hubble radius evolves as a 3 / 2 . 

We consider, for simplicity, a magnetic field whose typical correlation scale 
is as large as the Hubble radius at the electro-weak epoch when the temper- 
ature of the plasma was of the order of 100 GeV. This is roughly the regime 
contemplated by the considerations presented around Eq. (|4TJ)) . If the corre- 
lation scale evolves as the scale factor, the Hubble radius at the electroweak 
epoch (roughly 3 cm) projects today over a scale of the order of the astro- 
nomical unit. If inverse cascades are invoked, the correlation scale may grow, 
depending on the specific features of the cascade, up to 100 A.U. or even up to 
100 pc. In both cases the final scale is too small if compared with the typical 
scale of the gravitational collapse of the proto-galaxy. In Fig. [1] a particular 
model for the evolution of the correlation scale A (a) has been reported 0. 



Notice, as it will be discussed later, that the inverse cascade lasts, in principle, 
only down to the time of e + — e~ annihilation (see also thin dashed horizontal 
line in Fig. [1} since for temperatures smaller than T e +_ e - the Reynolds number 
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2.6 Inflationary magnetogenesis 

If magnetogenesis takes place inside the Hubble radius the main problem is 
therefore the correlation scale of the obtained seed field. The cure for this 
problem is to look for a mechanism producing magnetic fields that are coher- 
ent over large-scales (i.e. Mpc and, in principle, even larger). This possibility 
may arise in the context of inflationary models. Inflationary models may be 
conventional (i.e. based on a quasi-de Sitter stage of expansion) or unconven- 
tional (i.e. not based on a quasi-de Sitter stage of expansion) . Unconventional 
inflationary models are, for instance, pre-big bang models that will be dis- 
cussed in more depth in Section [3] 

The rationale for the previous statement is that, in inflationary models, the 
zero-point (vacuum) fluctuations of fields of various spin are amplified. Typi- 
cally fluctuations of spin and spin 2 fields. The spin 1 fields enjoy however 
of a property, called Weyl invariance, that seems to forbid the amplification 
of these fields. While Weyl invariance and its possible breaking will be the 
specific subject of the following subsection, it is useful for the moment to look 
at the kinematical properties by assuming that, indeed, also spin 1 field can 
be amplified. 

Since during inflation the Hubble radius is roughly constant (see Fig. [2]), 
the correlation scale evolves much faster than the Hubble radius itself and, 
therefore, large scale magnetic domains can naturally be obtained. Notice that, 
in Fig. [5] the (vertical) dashed lines illustrate the evolution of the Hubble ra- 
dius (that is roughly constant during inflation) while the full line denotes the 
evolution of the correlation scale. Furthermore, the horizontal (dashed) lines 
mark, from top to bottom, the end of the inflationary phase and the onset 
of the matter-dominated phase. This phenomenon can be understood as the 
gauge counterpart of the super-adiabatic amplification of the scalar and ten- 
sor modes of the geometry. The main problem, in such a framework, is to get 
large amplitudes for scale of the order of the Mpc at the onset of gravitational 
collapse. Models where the gauge couplings are effectively dynamical (break- 
ing, consequently, the Weyl invariance of the evolution equations of Abelian 
gauge modes) may provide rather intense magnetic fields. 

The two extreme possibilities mentioned above may be sometimes com- 
bined. For instance, it can happen that magnetic fields are produced by super- 
adiabatic amplification of vacuum fluctuations during an inflationary stage of 
expansion. After exiting the horizon, the gauge modes will reenter at different 
moments all along the radiation and matter dominated epochs. The spectrum 
of the primordial gauge fields after reentry will not only be determined by the 
amplification mechanism but also on the plasma effects. As soon as the mag- 
netic inhomogeneities reenter, some other physical process, taking place inside 
the Hubble radius, may be triggered by the presence of large scale magnetic 



drops below 1. This is the result of the sudden drop in the number of charged 
particles that leads to a rather long mean free path for the photons. 
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fields. An example, in this context, is the production of topologically non- 
trivial configurations of the hypercharge field (hypermagnetic knots) from a 
stochastic background of hypercharge fields with vanishing helicity 54, 531155] 
(see also [Si [53 [Si EH! EH)- 



SUPERADIABATIC AMPLIFICATION 




ICf'cm 



Fig. 2. Evolution of the correlation scale if magnetic fields would be produced by 
superadiabatic amplification during a conventional inflationary phase. The dashed 
vertical lines denote, in the present figure, the evolution of the Hubble radius while 
the full line denotes the evolution of the correlation scale (typically selected to 
smaller than the Hubble radius during inflation). 



2.7 Breaking of conformal invariance 

Consider the action for an Abelian gauge field in four-dimensional curved 
space-time 

1 ru 



S™ = — J d^xV-GF^F^. (48) 

Suppose, also, that the geometry is characterized by a conformally flat line 
element of Friedmann-Robertson- Walker type as the one introduced in Eq. 
([46]) . The equations of motion derived from Eq. ([48]) can be written as 



d„ V^GF^ = 0. (49) 



Using Eq. ([4*6"]) and recalling that v 7 — G = a 4 (r), we will have 

V=SF^ = a\r)^^ F) F afi =F^ (50) 

where the second equality follows from the explicit form of the metric. Equa- 
tion (|50p shows that the evolution equations of Abelian gauge fields are the 
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same in flat space-time and in a conformally flat FRW space-time. This prop- 
erty is correctly called Weyl invariance or, more ambiguously, conformal invari- 
ance. Weyl invariance is realized also in the case of chiral (massless) fermions 
always in the case of conformally flat space-times. 

One of the reasons of the success of inflationary models in making predic- 
tions is deeply related with the lack of conformal invariance of the evolution 
equations of the fluctuations of the geometry. In particular it can be shown 
that the tensor modes of the geometry (spin 2) as well as the scalar modes 
(spin 0) obey evolution equations that are not conformally invariant. This 
means that these modes of the geometry can be amplified and eventually af- 
fect, for instance, the temperature autocorrelations as well as the polarization 
power spectra in the microwave sky. 

To amplify large-scale magnetic fields, therefore, we would like to break 
conformal invariance. Before considering this possibility, let us discuss an even 
more conservative approach consisting in studying the evolution of Abelian 
gauge fields coupled to another field whose evolution is not Weyl invariant. 
An elegant way to achieve this goal is to couple the action of the hypercharge 
field to the one of a complex scalar field (the Higgs field) . The Abelian-Higgs 
model, therefore, leads to the following action 



where 2? M = <9 M — ieA^ and T^v = dt^A^ . Using Eq. (|16"]) into Eq. ([5"T]) and 
assuming that the complex scalar field (as well as the gauge fields) are not a 
source of the background geometry, the canonical action for the normal modes 
of the system can be written as 



where <P = a<j>\ = <9 M — ieA^ and F^ v = A u ] . From Eq. ([52]) it is clear 
that also when the Higgs field is massless the coupling to the geometry breaks 
explicitly Weyl invariance. Therefore, current density and charge density fluc- 
tuations will be induced. Then, by employing a Vlasov-Landau description 
similar the resulting magnetic field will be of the order of Bdoc ~ 10 ~ 4a 1% lc 
[52"] which is, by far, too small to seed any observable field even assuming, 
optimistically, perfect flux freezing and maximal efficiency for the dynamo ac- 
tion. The results of 62J disproved earlier claims (see [63] for a critical review) 
neglecting the role of the conductivity in the evolution of large-scale magnetic 
fields after inflation. 

The first attempts to analyze the Abelian-Higgs model in De Sitter space 
have been made by Turner and Widrow [66] who just listed such a possibil- 
ity as an open question. These two authors also analyzed different scenarios 
where conformal invariance for spin 1 fields could be broken in 4 space-time 
dimensions. Their first suggestion was that conformal invariance may be bro- 
ken, at an effective level, through the coupling of photons to the geometry 




(51) 




(52) 
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[69] . Typically, the breaking of conformal invariance occurs through products 
of gauge-field strengths and curvature tensors, i.e. 

-Lf^F^R^, J-R F^F» a g aP , -L Faf iF a PR (53) 
m z m z m -2 

where m is the appropriate mass scale; and R^ are the Riemann and 

Ricci tensors and R is the Ricci scalar. If the evolution of gauge fields is 
studied during phase of de Sitter (or quasi-de Sittter) expansion, then the 
amplification of the vacuum fluctuations induced by the couplings listed in 
Eq. ([53| is minute. The price in order to get large amplification should be, 
according to [66j , an explicit breaking of gauge-invariance by direct coupling 
of the vector potential to the Ricci tensor or to the Ricci scalar, i.e. 

RA^A", R^A^A". (54) 

In [66j two other different models were proposed (but not scrutinized in detail) 
namely scalar electrodynamics and the axionic coupling to the Abelian field 
strength. 

Dolgov [68j considered the possible breaking of conformal invariance due to 
the trace anomaly. The idea is that the conformal invariance of gauge fields is 
broken by the triangle diagram where two photons in the external lines couple 
to the graviton through a loop of fermions. The local contribution to the 
effective action leads to the vertex (iJ—g) 1+e F a pF al3 where e is a numerical 
coefficient depending upon the number of scalars and fermions present in the 
theory. The evolution equation for the gauge fields, can be written, in Fourier 
space, as 



A'l + ^HA' k + k 2 A k = 0, (55) 
8 

and it can be shown that only if e > the gauge fields are amplified. Further- 
more, only is e ~ 8 substantial amplification of gauge fields is possible. 

In a series of papers [701 [7TJ [72] the possible effect of the axionic coupling 
to the amplification of gauge fields has been investigated. The idea is here that 
conformal invariance is broken through the explicit coupling of a pseudo-scalar 
field to the gauge field (see Section 5), i.e. 

" F aP F^, (56) 



° m 87rM 

where F a/3 is the dual field strength and where numerical factor 

of order one. Consider now the case of a standard pseudoscalar potential, 
for instance m 2 ip 2 , evolving in a dc Sitter (or quasi-de Sitter space-time). It 
can be shown, rather generically, that the vertex given in Eq. (|56p leads to 
negligible amplification at large length-scales. The coupled system of evolution 
equations to be solved in order to get the amplified field is 

B" - V 2 B - ^V'V xB = 0, (57) 
ip" + 2Hip' + m 2 a 2 ip = 0, (58) 
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where B = a 2 B. From Eq. (|57[) . there is a maximally amplified physical 
frequency 



* w - M^ max= ^ m (59) 

where the second equality follows from ip ~ a~ 3 / 2 M cosmt (i.e. ?/>max ~ mM). 
The amplification for lu ~ w m ax is of the order of exp [ma e m/(27r.ff )] where 
H is the Hubble parameter during the de Sitter phase of expansion. From 
the above expressions one can argue that the modes which are substantially 
amplified are the ones for which w max 3> H. The modes interesting for the 
large-scale magnetic fields are the ones which are in the opposite range, i.e. 
w max "C H. Clearly, by lowering the curvature scale of the problem the pro- 
duced seeds may be larger and the conclusions much less pessimistic |72j . 

Another interesting idea pointed out by Ratra [67] is that the electro- 
magnetic field may be directly coupled to the inflaton field. In this case the 
coupling is specified through a parameter a, i.e. e av FapF" 13 where ip is the 
inflaton field in Planck units. In order to get sizable large-scale magnetic fields 
the effective gauge coupling must be larger than one during inflation (recall 
that cp is large, in Planck units, at the onset of inflation). 

In [73] it has been suggested that the evolution of the Abelian gauge cou- 
pling during inflation induce the growth of the two-point function of magnetic 
inhomogeneities. This model is different from the one previously discussed 
[67] . Here the dynamics of the gauge coupling is not related to the dynamics 
of the inflaton which is not coupled to the Abelian field strength. In particu- 
lar, rg(Mpc) can be as large as 10~ 12 . In [73] the MHD equations have been 
generalized to the case of evolving gauge coupling. Recently a scenario similar 
to [73] has been discussed in [74] . 

In the perspective of generating large scale magnetic fields Gasperini 
[75] suggested to consider the possible mixing between the photon and the 
graviphoton field appearing in supergravity theories (see also, in a related 
context [76]). The graviphoton is the massive vector component of the grav- 
itational supermultiplet and its interaction with the photon is specified by 
an interaction term of the type XF^G^" where G M „ is the filed strength of 
the massive vector. Large-scale magnetic fields with re(Mpc) > 10~ 34 can be 
obtained if A ~ 0(1) and for a mass of the vector m ~ 10 2 TeV. 

Bertolami and Mota [77] argue that if Lorentz invariance is spontaneously 
broken, then photons acquire naturally a coupling to the geometry which is 
not gauge-invariant and which is similar to the coupling considered in [66j . 



3 Why string cosmology? 

The moment has come to review my personal interaction with Gabriele 
Veneziano on the study of large-scale magnetic fields. While we had other 
15 joined papers with Gabriele (together with different combinations of au- 
thors) two of them [80] [81] (both in collaboration with Maurizio Gasperini) 
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are directly related to large-scale magnetic fields. Both papers reported in 
Refs. [SOI US] appeared in 1995 while I was completing my Phd at the theory 
division of CERN. 

My scientific exchange with Gabriele Veneziano started at least four years 
earlier and the first person mentioning Gabriele to me was Sergio Fubini. 
At that time Sergio was professor of Theoretical Physics at the University 
of Turin and I had the great opportunity of discussing physics with him at 
least twice a month. Sergio was rather intrigued by the possibility of getting 
precise measurements on macroscopic quantum phenomena like superfluidity, 
superconductivity, quantization of the resistivity in the (quantum) Hall effect. 
I started working, under the supervision of Maurizio Gasperini, on the spectral 
properties of relic gravitons and we bumped into the concept of squeezed state 



[82] . a generalization of the concept of coherent state (see, for instance, [85 ] l84 ] 

[55]). Sergio got very interested and, I think, he was independently thinking 
about possible applications of squeezed states to superconductivity, a topic 
that became later on the subject of a paper [86] . Sergio even suggested a review 
by Rodney Loudon [87] , an author that I knew already beacuse of his inspiring 
book on quantum optics [88]. Ref. [87] together with a physics report of B. L. 
Schumaker |89] was very useful for my understanding of the subject. Nowadays 
a very complete and thorough presentation of the intriguing problems arising 
in quantum optics can be found in the book of Leonard Mandel and Emil 



It is amusing to notice the following parallelism between quantum optics 
and the quantum treatment of gravitational fluctuations. While quantum op- 
tics deals with the coherence properties of systems of many photons, we deal, 
in cosmology, with the coherence properties of many gravitons (or phonons) 
excited during the time-evolution of the background fields. The background 
fields act, effectively, as a "pump field". This terminology, now generally ac- 
cepted, is exactly borrowed by quantum optics where the pump field is a laser. 
In the sixties and seventies the main problem of optics can be summarized by 
the following question: why is classical optics so precise? Put it into different 
words, it is known that the interference of the amplitudes of the radiation 
field (the so-called Young interferometry) can be successfully treated at a 
classical level. Quantum effects, in optics, arise not from the first -order inter- 
ference effects (Young interferometry) but from the second-order interference 
effects, i.e. the so-called Hanbury-Brown-Twiss interferometry [90 where the 
quantum nature of the radiation field is manifest since it leads, in the jargon 
introduced by Mandel [90] to light which is either bunched or anti-bunched. 
A similar problem also arises in the treatment of cosmological perturbations 
when we ask the question of the classical limit of a quantum mechanically 
generated fluctuation (for instance relic gravitons). 

The interaction with Sergio led, few years later, to a talk that I presented at 
the physics department of the University of Torino. The title was Correlation 
properties of many photons systems. I mentioned my interaction with Sergio 




Wolf [90]. 
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Fubini since it was Sergio who suggested that, eventually, I should talk to 
Gabriele about squeezed states. 

During the first few months of 1991, Gabriele submitted a seminal paper 
on the cosmological implications of the low-energy string effective action |91) . 
This paper, together with another one written in collaboration with Maur- 
izio Gasperini [92] represents the first formulation of pre-big bang models. A 
relatively recent introduction to pre-big bang models can be found in Ref. 
[93]. 

In [SHI [5T] it was argued that the string cosmological scenario provided 
by pre-big bang models [91] [92] would be ideal for the generation of large- 
scale magnetic fields. The rationale for this statement relics on two different 
observations: 

• in the low-energy string effective action gauge fields are coupled to the dila- 
ton whose expectation value, at the string energy scale, gives the unified 
value of the gauge and gravitational coupling; 

• from the mathematical analysis of the problem it is clear that to achieve 
a sizable amplification of large-scale magnetic fields it is necessary to have 
a pretty long phase where the gauge coupling is sharply growing in time 



Let us therefore elaborate on the two mentioned points. In the string 
frame the low-energy string effective action can be schematically written as 



In Eq. (jnOJ) the ellipses stand, respectively, for an expansion in powers of 
(A s /L) 2 and for an expansion in powers of the gauge coupling constant g 2 = 
e v . This action is written in the so-called string frame metric where the dilaton 
field if is coupled to the Einstein- Hilbert term. 

Concerning the action (|60|) few general comments are in order 

• the relation between the Planck and string scales depends on time and, 
in particular, £p = e v X 2 ; the present ratio between the Planck and string 
scales gives the value, i.e. g(r ) — e v °/ 2 — ^p(t )/A s ; 

• in four space-time dimensions the antisymmetric tensor field H^ va can be 
written in terms of a pseudo-scalar field, i.e. 



80] 



[Ml [SUES] 




+ 0(g 2 )+ ..(.60) 



\ivotp 




(61) 



In critical superstring theory the dilaton field must have a potential that van- 
ishes in the weak coupling limit (i.e. tp — > — oo). Moreover, from the direct 



32 Massimo Giovannini 



tests of Newton law at short distances it should also happen that the mass 
of the dilaton is such that m v > 10~ 4 . This requirement may be relaxed by 
envisaging non-perturbative mechanisms where the dilaton is effectively de- 
coupled from the matter fields and where a massless dilaton leads to observable 
violations of the equivalence principle. 

From the structure of the action (|60|) , Abelian gauge fields are amplified 
if the gauge coupling is dynamical. Consider, in fact, the equations of motion 
for the hypercharge field strength 



dJe-^V-GF^ J =0, (62) 

where F^ u = d^A v y In the Coulomb gauge where Ao — and V • A = 
the equation for the rescaled vector potential A^ = e^^A^ becomes, for each 
independent polarization and in Fourier space, 



A 



k z - 9 



A k = 0, (63) 



where, as usual, the prime denotes a derivation with respect to the conformal 
time coordinate. In Eq. (f6"3"]l k denotes the comoving wave-number From the 
structure of Eq. (|^|) there exist two different regimes. For k 2 ^> |<?(g _1 )"| the 
solution off Eq. (|63|) is oscillatory. In the opposite limit, i.e. k 2 <C ^(ff" 1 )"! 
the general solution can be written as 

Ci(A) c 2 (k) r , 

A k (T) = H — / g (t )dT , (64) 



9(r) g{r) 

where C\{k) and C2{k) are two arbitrary constants. These two constants can 
be fixed by imposing quantum mechanical initial conditions for r — > — oo. 
Thus, depending on the evolution of g(r) the Fourier amplitude Ak can be 
amplified. 

It can be shown [801 HI] that the the amplified magnetic energy density 
depends on the ratio between the value of the gauge coupling at the reentry 
and at the exit of the typical scale of the gravitational collapse, i.e. 

J_ dp B ^ fc 4 / £) 
p 7 d In k a 4 p 7 \ g, 



r(*) = fS-^r(P) ■ (65) 



The parameter r(k) measures the relative weight of the magnetic energy den- 
sity in units of the radiation background. To turn on the galactic dynamo in 
its simplest realization one should require that r(kc) > 10 -34 for a typical 
comoving wave-number corresponding to the typical scale of the gravitational 
collapse of the protogalaxy. As explained before, this requirement seems to be 
too optimistic in the light of the most recent understanding of the dynamo 
theory. The limit r(kc) > I0~ 24 seems more reasonable. 
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The fact that the gauge coupling must be sharply growing in order to 
produce large-scale magnetic fields fits extremely well with the pre-big bang 
dynamics where, indeed, the gauge coupling is expected to grow. The second 
requirement to obtain a phenomenologically viable mechanism for the ampli- 
fication of large-scale gauge fields turned out to be the existence of a pretty 
long stringy phase. 

The " stringy" phase is simply the epoch where quadratic curvature correc- 
tions start being important and lead to an effective dynamics where the dilaton 
field is linearly growing in the cosmic time coordinate (see |93j and references 
therein). Towards the end of the stringy phase the dilaton freezes to its (con- 
stant) value and the Universe gets dominated by radiation. One possibility for 
achieving the transition to radiation is represented by the back-reaction ef- 
fects of the produced particles |102j . In particular, the short wavelength modes 
play, in this context a crucial role. It is interesting that while the magnetic 
energy spectrum produced during the stringy phase is quasi-flat and the value 
of r(fcc) can be as large as 1CU 8 implying a protogalactic magnetic field of 
the order of 1CU 10 G. Under these conditions the dynamo mechanism would 
even be superfluous since the compressional amplification alone can amplify 
the seed field to its observed value. 

The results reported above may be "tested" in a framework where the pre- 
big bang dynamics is solvable. Consider, in particular, the situation where the 
evolution of the dilaton field as well as the one of the geometry is treated in 
the presence of a non-local dilaton potential [97 t M j l99 l HOO l [TUT] . 

In the Einstein frame description, the asymptotics of the (four-dimensional) 
pre-big bang dynamics can be written as |102j 



a(r) ~ a_4 / — - — > a_ = e 



2r ' V V3 



r 



W- = i, 9'- = -— , (66) 

zr t 



for t — > — co, and 



( ^ W2 / 2 (^3-l) 

a(r) " a+ Vw a+= 

, + ^o-ln2-^ln(^ f i) + V31n(^-), 

n + = ± ^ = ^, (67) 

for t — > +co. In Eqs. ((66)) and (|67p . H = a 1 /a and, as usual, the prime 
denotes a derivation with respect to r.The branch of the solution denoted by 
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minus describes, in the Einstein frame, an accelerated contraction, since the 
first derivative of the scale factor is negative while the second is positive. The 
branch of the solution denoted with plus describes, in the Einstein frame, a 
decelerated expansion, since the first derivative of the scale factor is positive 
while the derivative is negative. In both branches the dilaton grows and its 
derivative is always positive-definite (i.e. (p'± > ) as required by the present 
approach to bouncing solutions. The numerical solution corresponding to the 
asymptotics given in Eqs. ()66j) and (|67p is reported in Fig. [3] 




-100 -50 50 100 150 200 250 300 

X 

Fig. 3. The evolution of the scale factor in conformal time for a bouncing model 
regularized via non-local dilaton potential in the Einstein frame. 



In the Schrodinger description the vacuum state evolves, unitarily, to a 
multimode squeezed state, in full analogy with what happens in the case of 
relic gravitons [1031 11041 1105] . In the following the same process will be dis- 
cussed within the Heisenberg representation. The two physical polarizations 
of the photon can then be quantized according to the standard rules of quan- 
tization in the radiation gauge in curved space-times: 



a a / \ — ik-x i j-t oc a /_\* ik-; 



(68) 



and 



/d k 
j—^ a k , a e?n k (T)e-> k -*+al a e?II k (Tye 

OL V / - 

where ef(k) describe the polarizations of the photon and 

n k ( T ) = A' k (r) , [a k , a , aJ >/3 ] = S a0 5 ( ^ (k-p). 



* Ah-x 



(69) 



(70) 



The evolution equation for the mode functions will then be, in Fourier space, 
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A k = 0, 



(71) 



i.e. exactly the same equation obtained in (|63|) . The pump field can also be 
expressed as: 



g(g~ 



2 



The maximally amplified modes are then the ones for which 



(72) 
(73) 



The Fourier modes appearing in Eq. (|71[) have to be normalized while they are 
inside the horizon for large and negative r. In this limit the initial conditions 
provided by quantum mechanics are 



A(r) 



1 



n k {r) = -e 



(74) 



In the limit r — > +oo the positive and negative frequency modes will be 
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Fig. 4. The evolution of the mixing coefficients for k ~ fckmax in units of to. 



mixed, so that the solution will be represented in the plane wave orthonormal 
basis as 



A k (r) 
Mr) 



1 



~ 2 



c+{k)e- lkT +c_(fc)e lfcT 
c+{k)e- tkT -C-{k)e tkT 



(75) 



where c±(k) are the (constant) mixing coefficients. The following two relations 
fully determine the square modulus of each of the two mixing coefficients in 
terms of the complex wave- functions obeying Eq. (|7ip : 
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| c+ (fc)| 2 - Mfc)| 2 = i(A%n k - A k n* k ), (76) 

|c+(fc)| 2 + |c_(fc)| 2 - ^ (V fe | 2 + fc 2 |A| 2 ) ■ (77) 

After having numerically computed the time evolution of the properly nor- 
malized mode functions, Eqs. ([75]) and ([77)) can be used to infer the value 
of the relevant mixing coefficient (i.e. c_(fc)). Equation (|76[) is, in fact, the 
Wronskian of the solutions. If the second-order differential equation is writ- 
ten in the form (|7T|) . the Wronskian is always conserved throughout the time 
evolution of the system. Since, from Eq. (JTJJ, the Wronskian is equal to 1 
initially, it will be equal to 1 all along the time evolution. Thus, from Eq. ([75)) 
|c+(fc)| 2 = |c_(fc)| 2 + 1. The fact that the Wronskian must always be equal to 
1 is the measure of the precision of the algorithm. 

In Figs. [4] and [5] the numerical calculation of the spectrum is illustrated for 
different values of k. In Fig. 0the mixing coefficients are reported for modes 
k <C &max- In Fig. d] the mixing coefficients are reported for modes around 
t mH . Clearly, from Fig. [5] a smaller k leads to a larger mixing coefficient 
which means that the spectrum is rather blue. Furthermore by comparing 
the amplification of different modes it is easy to infer that the scaling law 
is |c+(fc)| 2 + |c_(fc)| 2 cx (fc/fc max )~™ 9 , with n g ~ 3.46, which is in excellent 
agreement with the analytical determination of the mixing coefficients leading 
to n g = 2V3 ~ 3.46[see below, Eq. (1551) ]. 

The second piece information that can be drawn from Fig. 3] concerns 
fcmax, whose specific value 

V5-0.5 

TO 

can be determined numerically for different values of to. 

For the value of £; max reported in Eq. (|75|) . the obtained mixing coefficient 
is 1, i.e. |c_(fc max )| ~ 1. According to Fig. [4] as we move from fc max to larger 
k, (|c+(fc)| 2 + |c_(fc)| 2 ) ~ (|c+(fc)| 2 - Mfc)| 2 ) implying that |c_(fc)| ~ 0. 
Moreover, from the left plot of Fig. [5] it can be appreciated that 

|c_(/c max )| 2 = 1, log(|c+(fc max )| 2 + Mfc max )| 2 ) = log3 ~ 0.477. (79) 

Thus the absolute normalization and slope of the relevant mixing coefficient 
can be numerically determined to be 

k ^ 2V ~ 3 



MAOIST— ■ ( 8 °) 



It can be concluded that Eq. ([80| is rather accurate as far as both the 
slope and the absolute normalization are concerned. The numerical estimates 
presented so far can be also corroborated by the usual analytical treatment 
based on the matching of the solutions for the mode functions before and 
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Fig. 5. The numerical estimate of the mixing coefficients in the case fcro <C 1. 



after the bounce. The evolution of the modes described by Eq. (|TTj) can be 
approximately determined from the exact asymptotic solutions given in Eqs. 
(|66|) and (|67|) . and implying that tp' ± ~ ±V3/t. Thus the solutions of Eq. ([7T]) 
can be obtained in the two asymptotic regimes, i.e. for r < — n 



and for t > t± 
A k ,+ (v) - 



if (n+l/2) 



e'2 



( 1 2) 

where i7„ (2) are Hankel functions of first and second kind whose related 
indices are 

v= — ^ — , M= — 2 — ' ( ' 

The time scale t± defines the width of the bounce and, typically t\ ~ tq. 
The phases appearing in Eqs. (f8Tj) and (|82|) are carefully chosen so that 



(81) 



r > -n, 
(82) 



lim „4j 



1 



-ikr 



2k 



(84) 



Using then the appropriate matching conditions 

A k ,-(-n) =A kt+ {n), 
A' k> _(-n)=A' kt+ (n) ) 

and defining x\ — fen, the obtained mixing coefficients are 



(85) 



c+(fe) 



1— x\e 
4 



i7r(y+At+l)/2 
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4 x\ p ■ 

+H<?) ( Xl )Hi% ( Xl )+ H%{ Xl )H^{ Xl ) 



(86) 



(87) 



satisfying the exact Wronskian normalization condition |c+(fc)| 2 — |c_(/c)| 2 = 
1. In the small argument limit, i.e. kr± ~ kro <C 1 the leading term in Eq. 
([57] leads to 

c_(fc) ~ — e ix ^-^ 2 Xi' , - v (y + n- l)r(u)r(u) (88) 

47T 

If we now insert the values given in Eq. (|83[) it turns out that c_(fc) ~ 
0.41 |fcTi|~^. The spectral slope agrees with the numerical estimate, as al- 
ready stressed. The absolute normalization cannot be determined from Eq. 
(|88p. where the small argument limit has already been taken. In order to deter- 
mine the absolute normalization the specific value of fc max ri has to be inserted 
in Eq. ([57)) . The result of this procedure, taking n ~ r is |c_(fc max )| 2 = 0.14, 
which is roughly a factor of 10 smaller than the interpolating formula given 
in Eq. (JgD}. 

The observation that a dynamical gauge coupling implies a viable mech- 
anism for the production of large-scale magnetic fields can be interesting in 
general terms and, more specifically, in the context of the pre-big bang models. 
In fact, in pre-big bang models, not only the fluctuations of the hypercharge 
field are amplified. In the minimal case we will have to deal with the fluctua- 
tions of the tensor [82, 106J and scalar |107j modes of the geometry and with 
the fluctuations of the antisymmetric tensor field [108, 109 . 

The amplified tensor modes of the geometry lead to a stochastic back- 
ground of gravitational waves (GW) with violet spectrum both in the GW 
amplitude and energy density. In Fig. [6] the GW signal is parametrized in 
terms of the logarithm of J7qw = pGw/pc, i-e. the fraction of critical energy 
density present (today) in GW. On the horizontal axis of Fig.[6]the logarithm 
of the present (physical) frequency v is reported. In conventional inflation- 
ary models, for v > 10~ 16 Hz, 4?gw, is constant (or slightly decreasing) as 
a function of the present frequency. In the case of string cosmological mod- 
els f^GW °c is 3 Inv, which also implies a steeply increasing power spectrum. 
This possibility spurred various experimental groups to analyse possible di- 
rects limits on the scenario arising from specific instruments such as resonant 
mass detectors |110j and microwave cavities |111[ 1112] . These attempts are 
justified since the signal of pre-big bang models may be rather strong at high 
frequencies and, anyway, much stronger than the conventional inflationary 
prediction 

The sensitivity of a pair of VIRGO detectors to string cosmological gravi- 
tons has been specifically analysed |113j with the conclusion that a VIRGO 
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Fig. 6. The spectrum of relic gravitons from various cosmological models presented 
in terms of /i 2 j9gw ■ 

pair, in its upgraded stage, can certainly probe wide regions of the parameter 
space of these models. If we maximize the overlap between the two detec- 
tors [113j or if we reduce (selectively) the pendulum and pendulum's internal 
modes contribution to the thermal noise of the instruments, the visible region 
(after one year of observation and with SNR = 1) of the parameter space will 
get even larger. Unfortunately, as in the case of the advanced LIGO detectors, 
the sensitivity to a flat J?gw will be irrelevant for ordinary inflationary models 
also with the advanced VIRGO detector. It is worth mentioning that growing 
energy spectra of relic gravitons can also arise in the context of quintessential 
inflationary models [1 141 1115] . In this case J^gw v In 2 v (see [115] for a full 
discussion) . 

The spectra of gravitational waves have features that are, in some sense, 
complementary to the ones of the large-scale magnetic fields. The parameter 
space leading to a possible signal of relic (pre-big bang) gravitons with wide- 
band interferometers has only a small overlap with the region of the parameter 
space leading to sizable large-scale magnetic fields. This conclusion can be 
evaded if the coupling of the dilaton to the hypercharge field is, in the action, 
of the type e - ^ F flu F lil/ [116] where the parameter (3 has values 1 and 1/2, 
respectively, for heterotic and type I superstrings. In particular, in the case 
(3 = 1/2, it is possible to find regions where both large-scale magnetic fields 
and relic gravitons are copiously produced. 

Let us finally discuss the scalar fluctuations of the geometry. The spec- 
trum of the scalar modes is determined by the spectrum of the Kalb-Ramond 
axion(s). If the axions would be neglected, the spectrum of the curvature fluc- 
tuations would be sharply increasing, or as we say in the jargon, the spectrum 
would be violet in full analogy with the spectrum of the tensor modes of the 
geometry. This result |107j has been recently analyzed in the light of a recent 
controversy (see [57J|5H]) and references therein). 
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If the Kalb-Ramond axions are consistently included in the calculation, 
it is found that the large-scale spectrum of curvature perturbations becomes 
flat [109j and essentially inherits the spectrum of the Kalb-Ramond axions. 
If the axions decay (after a phase of coherent oscillations) the curvature per- 
turbations will be adiabatic as in the case of conventional inflationary models 
but with some important quantitative differences 109J since, in this case, the 
CMB normalization is explained in terms of the present value of the string 
curvature scale and in terms of the primordial slope of the axion spectrum. 

4 Primordial or not primordial, this is the question... 

While diverse theoretical models for the origin of large-scale magnetism can 
certainly be questioned on the basis of purely theoretical considerations, di- 
rect observations can tell us something more specific concerning the epoch 
of formation of large-scale magnetic fields. It would be potentially useful to 
give some elements of response to the following burning question: are really 
magnetic fields primordial? 

The plan of the present section is the following. In Subsect. 14.11 differ- 
ent meanings of the term primordial will be discussed. It will be argued that 
CMB physics can be used to constrain large-scale magnetic fields possibly 
present prior to matter-radiation equality. In Subsect. 14.21 the scalar CMB 
anisotropies will be specifically discussed by deriving the appropriate set of 
evolution equations accounting for the presence of a fully inhomogeneous mag- 
netic field. In Subsect. 14. 3l the evolution of the different species composing the 
pre-decoupling plasma will be solved, in the tight-coupling approximation and 
in the presence of a fully inhomogeneous magnetic field. Finally Subsect. 14.41 
contains various numerical results and a strategy for parameter extraction. 

4.1 Pre-equality magnetic fields 

The term primordial seems to have slightly different meanings depending on 
the perspective of the various communities converging on the study of large- 
scale magnetic fields. Radio-astronomers have the hope that by scrutinizing 
the structure of magnetic fields in distant galaxies it would be possible, in the 
future, to understand if the observed magnetic fields are the consequence of a 
strong dynamo action or if their existence precedes the formation of galaxies. 

If the magnetic field does not flips its sign from one spiral arm to the 
other, then a strong dynamo action can be suspected |117j . In the opposite 
case the magnetic field of galaxies should be primordial i.e. present already at 
the onset of gravitational collapse. In this context, primordial simply means 
protogalactic. An excellent review on the evidence of magnetism in nearby 
galaxies can be found in [118j . In Fig. [7] a schematic view of the Milky Way is 
presented. The magnetic field follows the spiral arm. There have been claims, 
in the literature, of 3 to 5 field reversals. The arrows in Fig. [7] indicate one of 
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Fig. 7. The schematic map of the MW is illustrated. Following |119| the origin of 
the two-dimensional coordinate system are in the Galactic center. The two large 
arrows indicate one of the possible (3 or 5) field reversals observed so far. 



the possible field reversals. One reversal is certain beyond any doubt. Another 
indication that would support the primordial nature of the magnetic field of 
galaxies would be, for instance, the evidence that not only spirals but also 
elliptical galaxies are magnetized (even if the magnetic field seems to have 
correlation scale shorter than in the case of spirals). Since elliptical galaxies 
have a much less efficient rotation, it seems difficult to postulate a strong 
dynamo action. We will not pursue here the path of specific astrophysical 
signatures of a truly pre-galactic magnetic field and e refer the interested 
reader to [TT7lflT8] . 

As a side remark, it should also be mentioned that magnetic fields may play 
a role in the analysis of rotation curves of spiral galaxies. This aspect has been 
investigated in great depth by E. Battaner, E. Florido and collaborators also 
in connection with possible effects of large-scale magnetic fields on structure 
formation [1201 USD DM [123] (see also [124] and references therein). 

The large-scale magnetic fields produced via the parametric amplification 
of quantum fluctuations discussed earlier in the present lecture may also be 
defined primordial but, in this case, the term primordial has a much broader 
signification embracing the whole epoch that precedes the equality between 
matter and radiation taking place, approximately, at a redshift z cq = 3230 for 
h 2 fl m Q = 0.134 and h 2 n i0 — 4.15 x 10 -5 . Consequently, large-scale magnetic 
fields may affect, potentially, CMB anisotropies [19] . Through the years, vari- 
ous studies have been devoted to the effect of large-scale magnetic fields on the 
vector and tensor CMB anisotropies [1281 1129j (see also [125] and references 
therein for some recent review articles). 

The implications of fully inhomogeneous magnetic fields on the scalar 
modes of the geometry remain comparatively less explored. By fully inho- 
mogeneous we mean stochastically distributed fields that do not break the 
spatial isotropy of the background |1251 1126j . 
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CMB anisotropies are customarily described in terms of a set of carefully 
chosen initial conditions for the evolution of the brightness perturbations of 
the radiation field. One set of initial conditions corresponds to a purely adi- 
abatic mode. There are, however, more complicated situations where, on top 
of the adiabatic mode there is also one (or more) non-adiabatic mode(s). A 
mode, in the present terminology, simply means a consistent solution of the 
governing equations of the metric and plasma fluctuations, i.e. a consistent 
solution of the perturbed Einstein equations and of the lower multipoles of 
the Boltzmann hierarchy 

The simplest set of initial conditions for CMB anisotropies, implies, in a 
/1CDM framework, that a nearly scale-invariant spectrum of adiabatic fluctu- 
ations is present after matter-radiation equality (but before decoupling) for 
typical wavelengths larger than the Hubble radius at the corresponding epoch 

urn 

It became relevant, through the years, to relax the assumption of exact 
adiabaticity and to scrutinize the implications of a more general mixture of 
adiabatic and non-adiabatic initial conditions (see [132[ I133[ I134j and refer- 
ences therein). In what follows it will be argued, along a similar perspective, 
that large-scale magnetic fields slightly modify the adiabatic paradigm so that 
their typical strengths may be constrained. To achieve such a goal, the first 
step is to solve the evolution equations of magnetized cosmological pertur- 
bations well before matter-radiation equality. The second step is to follow 
the solution through equality (and up to decoupling). On a more technical 
ground, the second step amounts to the calculation of the so-called trans- 
fer matrix [135j whose specific form is one of the the subjects of the present 
analysis. 

4.2 Basic Equations 

Consider then the system of cosmological perturbations of a flat Friedmann- 
Robertson- Walker (FRW) Universe, characterized by a conformal time scale 
factor a(r) (see Eq. (|46|) ). and consisting of a mixture of photons, baryons, 
CDM particles and massless neutrinos. In the following the basic set of equa- 
tions used in order to describe the magnetized curvature perturbations will be 
introduced and discussed. The perspective adopted here is closely related to 
the recent results obtained in Refs. [136, 137J (see also [138, 139 for interesting 
developments) . 

In the conformally Newtonian gauge |140[ I141[ I142|. 11431 H44j , the scalar 
fluctuations of the metric tensor G M „ = a?(r)r\ iiv are parametrized in terms 
of the two longitudinal fluctuations i.e. 

SGoo = 2a 2 (r)^{T : x), 5G l3 = 2a 2 (r)^(r, x)8 ih (89) 

where Sij is the Kroeneker S. While the spatial curvature will be assumed to 
vanish, it is straightforward to extend the present considerations to the case 
when the spatial curvature is not negligible. 
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In spite of the fact that the present discussion will be conducted within 
the conformally Newtonian gauge, it can be shown that gauge-invariant de- 
scriptions of the problem are possible [137] . Moreover, specific non-adiabatic 
modes (like the ones related to the neutrino system) may be more usefully 
described in different gauges (like the synchronous gauge). The rationale for 
the last statement is that the neutrino isocurvature modes may be singular 
in the conformally Newtonian gauge. These issues will not be addressed here 
but have been discussed in the existing literature (see, for instance, |143|I144| 
and references therein). Furthermore, for the benefit of the interested reader 
it is appropriate to mention that the relevant theoretical tools used in the 
present and in the following paragraphs follows the conventions of a recent 
review [144] . 

Hamiltonian and momentum constraints 

The Hamiltonian and momentum constraints, stemming from the (00) and 
(Oi) components of the perturbed Einstein equations are: 

VV - m{U<t> + ip') = 4irGa 2 [5 Pt + 6p B ], (90) 

V- (E x B)~ 



(Pt + Pt)0t + 



47T0 4 



(91) 



where TL — a' /a and the prime denotes a derivation with respect to the 
conformal time coordinate r. In writing Eqs. (|90|) and ([9l"|) the following set 
of conventions has been adopted 

Spt(r, x) = ^p 7 (r, x) + 5p v (r, x) + Sp c (r, x) + Sp h (r, x), (92) 

(Pt + pt)6 t (r,x) = (p 7 +p 7 )6» 7 (r, x) + (jp v + p u )6 v (r,x) 

+ (p c + Pc)0 c (t, x) + (pt + p h )6 h (T, x). (94) 



Concerning Eqs. (|92[) . ([93)1 and (|94|) the following comments are in order: 

• in Eq. (|92[) the total density fluctuation of the plasma, i.e. Spt(r, x) receives 
contributions from all the species of the plasma; 

• in Eq. (|93[) the fluctuation of the magnetic energy density Spb(t, x) is 
quadratic in the magnetic field intensity; 

• in Eq. (f9"4")l 9 t (T,x) — div\ is the divergence of the total peculiar velocity 
while 8 7 (t,x), 9 v (t, x), 9 c (r,x) and 0b(T, x) are the divergences of the 
peculiar velocities of each individual species, i.e. photons, neutrinos, CDM 
particles and baryons. 

The second term appearing at the right hand side of Eq. (f9"Tj) is the diver- 
gence of the Poynting vector. In MHD the Ohmic electric field is subleading 
and, in particular, from the MHD expression of the Ohm law we will have 
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Ex B ~ 



(V xB)xB 



(95) 



Aire 



Since the Universe, prior to decoupling, is a very good conductor, the ideal 
MHD limit can be safely adopted in the first approximation (see also [145] ); 
thus for (t — > oo (i.e. infinite conductivity limit) the contribution of the Poynt- 
ing vector vanishes. In any case, even if a would be finite but large, the second 
term at the right hand side of Eq. (|9"Tj) would be suppressed in comparison 
with the contribution of the divergence of the total velocity field. 

The total (unperturbed) energy density and pressure of the mixture, i.e. 



determine the evolution of the background geometry according to Friedmann 
equations: 



Notice that in Eq. (|96| the contribution of the cosmological constant has 
been included. If the dark energy is parametrized in terms of a cosmological 
constant (i.e. pa = —Pa), then, 5p' A = 0. Furthermore, the contribution of pa 
to the background evolution is negligible prior to decoupling. Slightly different 
situations (not contemplated by the present analysis) may arise if the dark 
energy is parametrized in terms of one (or more) scalar degrees of freedom 
with suitable potentials. 

Dynamical equation and anisotropic stress (es) 

The spatial components of the perturbed Einstein equations, imply, instead 



Equation (| 1 00[) contains, as source terms, not only the total fluctuation of the 
pressure of the palsma, i.e. Spt, but also 



Pt = Pi + Pv + Pc + Ph + Pa, 

Pt = P-y + Pv + Pc + Pb + PA- 



(96) 




(97) 

(98) 
(99) 



H 2 -W = 4nGa 2 {p t + Pt ), 
p[ + 3W(p t +Pt) = 0. 



ijy + n(<p' + 2y/) + {h 2 + 2H')<p + -v 2 (0 - y) 5{ 

-icW^-VO =^Ga 2 (S Pt + 6 PB )Sf - nf - fl{ . (100) 



Spb{t, x) 



B 2 {x) 



Sp B (T,x) 
3 



(101) 



24tto 4 (t) 




(102) 
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Moreover, in Eq. (|100|) , II?(t,x) is the anisotropic stress of the fluid. As it 
will be mentioned in a moment (and later on heavily used) the main source of 
anisotropic stress of the fluid is provided by neutrinos which free-stream from 
temperature smaller than the MeV. Notice that both the anisotropic stress of 
the fluid, i.e. W(t,x) and the magnetic anisotropic stress, i.e. 77f (r, x), are, 
by definition, traceless. 

Using this last observation, Eq. (|100p can be separated into two indepen- 
dent equations. Taking the trace of Eq. (|100|) we do get 

<tp" + H{4>' + 2ip') + (2H' + n 2 )4> + \v 2 {4> -ip) = 4irGa 2 {5 Pt + 6 Pb ). (103) 

By taking the difference between Eq. (| 1 00|) and Eq. (| 1 03|) the following (trace- 
less) relation can be obtained: 

cW(0 - V) - \siV 2 (4> -il>) = 8irGa 2 {nj + 77/). (104) 

By applying the differential operator djd 1 to both sides of Eq. (|104j) we do 
obtain the following interesting relation: 

V 4 (0 - V) = l2vGa 2 [{p v + f>u)V 2 a v + (p 7 + p 7 )V 2 cr B ], (105) 

where the parametrization 

d^nl = (p v + p„)V 2 a u , = (p 7 + i o 7 )V 2 a B , (106) 

has been adopted. In Eq. (|105[) <j v {t, x) is related with the quadrupole moment 
of the (perturbed) neutrino phase-space distribution. In Eq. (|105p <7b(t, x ) 
parametrizes the (normalized) magnetic anisotropic stress. It is relevant to 
remark at this point that in the MHD approximation adopted here the two 
main sources of scalar anisotropy associated with magnetic fields can be 
parametrized in terms of (Jb{t, x) and in terms of the dimensionless ratio 

O b (t,x) = 6 J^. (107) 

Since both J7b(t, x) and ub(t, x) are quadratic in the magnetic field intensity 
a non-Gaussian contribution may be expected. f2 B {r, x) is the magnetic energy 
density referred to the photon energy density and it is constant to a very good 
approximation if magnetic flux is frozen into the plasma element. 

There is, in principle, a third contribution to the scalar problem coming 
from magnetic fields. Such a contribution arises in the evolution equation of 
the photon-baryon peculiar velocity and amounts to the divergence of the 
Lorentz force. While the mentioned equation will be derived later in this 
section, it is relevant to point out here that the MHD Lorentz force can be 
expressed solely in terms of <tb {t, x) and X?b (t, x) . In fact a well known vector 
identity stipulates that 
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d i B j d j B i = V • [(V x 73) x 73] + ^W 2 B 2 . (108) 
From the definition of ctb in terms of 77/, i.e. Eq. (|106|) . it is easy to show that 

vV b = r R 3 4 9 i B j d^B i - iv 2 ^ B . (109) 

Using then Eq. (TTU5|) into Eq. (TTUg|) and recalling that 

4ttV ■ [J x 73] = V • [(V x B) x B], (110) 

we obtain: 

V 2 ob = t^t^V • [(V x 73) x 73] + (111) 
Curvature perturbations 

Two important quantities must now be introduced. The first one, convention- 
ally denoted by £, is the density contrast on uniform curvature hypersurfaces 
B i.e. 

C= -^-tt (fo + frB) . (112) 
Pi 

The definition (| 1 1 2|) is invariant under infinitesimal coordinate transforma- 
tions. In fact, while SpB is automatically gauge-invariant (since the magnetic 
field vanishes at the level of the background) t/j and Spt transform as [144j 

t/j — > i/j = tp + He, 

Spt -» 6p t - p[e. (113) 

for 

r -> f = t + e° 

a; 4 -» £ 4 = a; 2 + d z e. (114) 
Recalling Eq. Eq. (|112[) can also be written as 



C = ^ + |^. (115) 
3(Pt +Pt) 

The second variable we want to introduce, conventionally denoted by 1Z is 
the curvature perturbation on comoving orthogonal hypersurfaces 

EE i.e. 



11 Since, as it will be discussed, Q is gauge-invariant, we can also interpret it as the 
curvature fluctuation on uniform density hypersurfaces, i.e. the fluctuation of the 
scalar curvature on the hypersurface where the total density is uniform. 

12 It is clear, from the definition (|116[1 that the second term at the right hand side 
is proportional, by the momentum constraint (|91[) , to the total peculiar velocity 
of the plasma which is vanishing on comoving (orthogonal) hypersurfaces. 
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^ = ' V " 7?-M' ■ (116) 

Inserting Eq. (|1 15[) and (|116[) into Eq. (|90[) . the Hamiltonian constraint takes 
then the form 

127rGo 2 (p t +fh) 

Equation l|117|) is rather interesting in its own right and it tells that, in the 
long wavelength limit, 

C~n + 0{k 2 T 2 ). (118) 

When the relevant wavelengths are larger than the Hubble radius (i.e. kr <C 1) 
the density constrast on uniform curvature hypersurfaces and the curvature 
fluctuations on comoving orthogonal hypersurfaces coincide. Since the ordi- 
nary Sachs- Wolfe contribution to the gauge-invariant temperature fluctuation 
is dominated by wavelengths that are larger than the Hubble radius after mat- 
ter radiation equality (but before radiation decoupling), the calculation of £ 
(or TV), in the long wavelength limit, will essentially give us the Sachs- Wolfe 
plateau. 

A remark on the definition given in Eq. (|112[) is in order. The variable 
£ must contain the total fluctuation of the energy density. This is crucial 
since the Hamiltonian constraint is sensitive to the total fluctuation of the 
energy density. If the magnetic energy density SpB is correctly included in the 
definition of £, then the Hamiltonian constraint (|117|) maintains its canonical 
form. 

Equations (fTT7|) and pT8|) can be used to derive the appropriate transfer 
matrices, allowing, in turn, the estimate of the Sachs- Wolfe plateau. For this 
purpose it is important to deduce the evolution equation for £. The evolution of 
£ can be obtained from the evolution equation of the total density fluctuation 
which reads, in the conformally Newtonian gauge, 

E J 

5p[ - 3V>'(Pt + pt) + (Pt + PtWt + 3H(5pt + 6pt) + 3H5 Pnad = — j-. (119) 
The technique is now rather simple. We can extract 5pt from Eq. ()1 15[) 

apt = 3(pt+Pt)(C + V0-*/>B- ( 12 °) 

Inserting Eq. Q120|) into Eq. (|119[) we get to the wanted evolution equation 
for £. Before doing that it is practical to discuss the case when the relativis- 
tic fluid receives contributions from different species that are simultaneously 
present. In the realistic case, considering that the cosmological constant does 
not fluctuate, we will have four different species. 

For deriving the evolution equation of £, it is practical (and, to some 
extent, conventional) to separate the pressure fluctuation into an adiabatic 
component supplemented by a non-adiabatic contribution: 
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to = (P) to + ) fc. (121) 



toy, v^/pt 

In a relativistic description of gravitational fluctuations, the pressure fluctu- 
ates both because the energy density fluctuates (first term at the right hand 
side of Eq. (|121[) ) of because the specific entropy of the plasma, i.e. q fluctu- 
ates (first term at the right hand side of Eq. (|121|l ). The subscripts appearing 
in the two terms at the righ-hand side of Eq. (|12ip simply mean that the two 
different variations must be taken, respectively, at constant <; (i.e. 5<; — 0) and 
at constant pt (i.e. dpt = 0). 

Here is an example of the usefulness of this decomposition. Consider, for 
instance, a mixture of CDM particles and radiation. In this case the coefficient 
of the first term at the right hand side of Eq. <| 1 2 1 j) can be written as 



Sp t \ if Sp_ 



to/c 3\Sp c + 5p 



(122) 



where we simply used the fact that Sp r — <5p r /3 and that 5pt = 5pr + Sp c . 
Now, the quantity appearing in Eq. (|122[) must be evaluated at constant ?, 
i.e. for 5? = 0. The specific entropy, in the CDM radiation system, is given 
by ? = T 3 /n c where T is the temperature and n c is the CDM concentration. 
The relative fluctuations of the specific entropy can then be defined and they 
are 

5= *£ = 3to_to 

C 4 Pr Pc 

where it has been used that p r ~ T 4 and that p c ~ m n c (m is here the 
typical mass of the CDM particle). Requiring now that S = we do get 5p c — 
(i/A)(p c /p 1 )8p r . Thus, inserting Sp c into Eq. (|122p . the following relation can 
be easily obtained: 

to/ ? 3(3p c + 4p r ) p[ 

The second and third equalities in Eq. (|124p follow from the definition 
of the total sound speed for the CDM-radiation system. This occurrence is 
general and it is not a peculiarity of the CDM-radiation system so that we 
can write, for an arbitrary mixture of relativistic fluids: 

to \ _ p' t _ 2 



toy\ Pt 

The definition of relative entropy fluctuation proposed in Eq. (|123p is invariant 
under infinitesimal gauge transformations [144] and it can be generalized by 
introducing two interesting variables namely 

Cr = -V-W% Cc = -V>-W% (126) 

Pr Pc 
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Using the continuity equations for the CDM and for radiation, i.e. p[ = — £Hp t 
and p' c = —3TCp c , Eq. (|126p can be also written as 

C = -^+|, Cc = -1> + j, (127) 

where 5 r = Sp r /p r and S c = Sp c /p c . Thus, using Eq. (|127p . the relative fluc- 
tuation in the specific entropy introduced in Eq. (| 1 23|) can also be written 
as 

S = -3(Cc - Cr)- (128) 

It is a simple exercise to verify that Eqs. (jl23[> and (|128|) have indeed the 
same physical content. 

Up to now the coefficient of the first term at the right-hand side of Eq. 
(I121|) has been computed. Let us now discuss the second term appearing at 
the right hand side of Eq. (|121jl . Conventionally, the whole second term is 
often denoted by Sp aa d, i.e. non-adiabatic pressure variation. From Eq. (1123|) 
defining the relative fluctuation in the specific entropy, i. e. S = 5?/?, the 
following equation can be written: 

Now, S must be evaluated, inside the round bracket, for Spt — 0. The result 
will be 

'Spt\ _4 p cPl 



S J pt 3 3,0c + 4p r 

Recalling the definition of sound speed and using Eq. (|130|) into Eq. (|129[) . we 
do get 

Spnad = C 2 s p c S. (131) 

If the mixture of fluids is more complicated the discussion presented so 
far can be easily generalized. If more than two fluids are present, we can still 
separate, formally, the pressure fluctuation as 

Spt = c 2 s 5p t + 6p nad . (132) 

However, if more than two fluids are present, the non-adiabatic pressure den- 
sity fluctuation has a more complicated form that reduces to the one previ- 
ously computed in the case of two fluids: 

6pnad = 6777 5- P ' { P ^ ~ C ' j)5ij ' 

* ij 

5 ij = -3(Ci-Cj), 4 = ^. (133) 

where <Si j are the relative fluctuations in the entropy density that can be 
computed in terms of the density contrasts of the individual fluids. The indices 
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i and j run over all the components of the plasma. Assuming a plasma formed 
by photons, neutrinos, baryons and CDM particles we will have that various 
entropy fluctuations are possible. For instance 

5 7C = -3(C 7 -Cc), Sy V = -3(C 7 - ■- ( 134 ) 

where the ellipses stand for all the other possible combinations. From the 
definition of relative entropy fluctuations it appears that Sy V = —S V1 . Finally, 
with obvious notations, while c 2 denotes the total sound speed, c 2 ; and c 2 ; 
denote the sound speeds of a generic pair of fluids contributing S-^ to 5p nat \, 
i.e. f 

c s = — 7' c si — — n c sj — ~7- (135) 
Pt Pi Pj 

In the light of Eq. (|134p . also the physical interpretation of Eq. (|132|) becomes 
more clear. The contribution of <5p na d arises because of the inherent multiplic- 
ity of fluid present in the plasma. Thanks to Eq. Q132|) using Eq. (|120|) in Eq. 
(|119|) we can obtain the evolution equation for £ which becomes 

H H ( 2 1 



C = —Spn* + —— K - - )Spb - -f (136) 

Pt+Pt Pt+Pt\ 3/ 3 

The evolution equation for 1Z can also be directly obtained by taking the 
first time derivative of Eq. (|117p . i.e. 

C yj | V ^ y | U{ic " + ^ (137) 

^ ^ 127rGa 2 (p t +Pt) UirGa? ( Pt + p t ) ' { ' 

By now inserting Eq. (|13T|) into Eq. (|136|) and by using the momentum con- 
straint of Eq. (j9T|) to eliminate 8t we do get the following expression: 

h n ( 2 i 



^' = ; dPnad H ; Us - o ) S P B 

Pt + Pt Pt+Pt\ 3/ 

+ m_^„ 2 /„ ; rr- ( ld8 -) 



47rGa 2 (p t + Pt) 127rGa 2 (p t + Pt) ' 

It could be finally remarked that Eq. (1138j) can be directly derived from Eq. 
(I103P . For this purpose The definition (I116P can be derived once with respect 
to t. The obtained result, once inserted back into Eq. (|103|) reproduces Eq. 

(HMD- 



4.3 Evolution of different species 

Up to now the global variables defining the evolution of the system have been 
discussed in a unified perspective. The evolution of the global variables is 
determined by the evolution of the density contrasts and peculiar velocities of 
the different species. Consequently, in the following paragraphs, the evolution 
of the different species will be addressed. 
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Photons and baryon 

The evolution equations of the lowest multipoles of the photon-baryon system 
amount, in principle, to the following two sets of equations: 



and 



8' b = 3V/ - 



H9 h = -V z 



V • [ J x B] A Pl 



a 4 ph 



+ o — an e x e CTT{9j - Ob), 

o Ph 



(139) 
(140) 



5; = 40 - ~0 7 , 
„, V 2 (L - 



an c x c a T ( 



(141) 
(142) 



Equation (| 140[) contains, as a source term, the divergence of the Lorentz force 
that can be expressed in terms of 0b(t, x) and f2-a{r,x), as already pointed 
out in Eqs. (TTTTj) . 

At early times photons and baryons are tightly coupled by Thompson 
scattering, as it is clear from Eqs. (|140p and (1142[) where ctt denotes the 
Thompson cross section and n c x c the concentration of ionized electrons. To 
cast light on the physical nature of the tight coupling approximation let us 
subtract Eqs. fH^]) and fliP]) . The result will be 



3 Ph 



V 2 cL 



V ■ [J x B] 



a 4 pb 



(143) 



From Eq. (| 143[) it is clear that any deviation of {0 1 — 9b) swiftly decays away. 
In fact, from Eq. (|143p . the characteristic time for the synchronization of the 
baryon and photon velocities is of the order of (:Eerc.eOT) -1 which is small 
compared with the expansion time. In the limit <tt — > 00 the tight coupling is 
exact and the photon-baryon velocity field is a unique physical entity which 
will be denoted by (9 7 b- From the structure of Eq. (| 143[) . the contribution of 
the magnetic fields in the MHD limit only enters through the Lorentz force 
while the damping term is always provided by Thompson scattering. 

To derive the evolution equations for the photon-baryon system in the 
tight coupling approximation we can add Eqs. (|140[) and (| 142|) taking into 
account that 9b — 9^ = # 7 b- Of course, also the evolution equations of the 
density contrasts will depend upon 9 7 b- Consequently the full set of tightly 
coupled evolution equations for the photon-baryon fluid can be written as: 



S' h = 3tp' - 6 7 b, 
, HRb 

^ b+ (TTi4) 



7 b 



7 7 b 



4(l + i? b ) 



3 V.[JXB] 



4aV 7 (l + i?b) 



(144) 
(145) 

(146) 
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where 

3 Pb (r) / 698 \(h*f2 b \ 

The set of equations (|144|) , (|145|) and (|146|) have to be used in order to obtain 
the correct initial conditions to be imposed on the evolution for the integration 
of the brightness perturbations. 

If we assume, effectively, that ctt — ► oo we are working to lowest order 
in the tight coupling approximation. This means that the CMB is effectively 
isotropic in the baryon rest frame. To discuss CMB polarization in the presence 
of magnetic helds one has to go to higher order in the tight coupling expansion. 
However, as far as the problem of initial conditions is concerned, the lowest 
order treatment suffices, as it will be apparent from the subsequent discussion. 



Neutrinos 

After neutrino decoupling the (perturbed) neutrino phase space distribution 
evolves according to the collisionless Boltzmann equation. This occurrence im- 
plies that to have a closed system of equations describing the initial conditions 
it is mandatory to improve the fluid description by adding to the evolution of 
the monopole (i.e. the neutrino density contrast) and of the dipole (i.e. the 
neutrino peculiar velocity) also the quadrupole, i.e. the quantity denoted by 
tj v and appearing in the expression of the anisotropic stress of the fluid (see 
Eqs. (HUSD and (TITO)) ). 

The derivation of the various multipoles of the perturbed neutrino phase 
space distribution is a straightforward (even if a bit lengthy) calculation and 
it has been performed, for the set of conventions employed in the present 
lecture, in Ref. [144] . The result is, in Fourier space, 

K = W - \e v , (148) 

0'„ = ^V 2 ^ + k 2 q) - k 2 a Ul (149) 

o' v = ^-6 v -^-kT vZ . (150) 
15 10 

In Eq. (|150[) T v z is the octupole of the (perturbed) neutrino phase space 
distribution. The precise relation of the multipole moments of T v with the 
density contrast and the other plasma quantities is as follows: 

5 V = Fvo, V = -kF v i, <T y = (151) 

For multipoles larger than the quadrupole, i.e. I > 2 the Boltzmann hierarchy 
reads: 

Ki = 2^X1^-1) - V + V^t+i)]- ( 152 ) 



Magnetic fields, strings and cosmology 



53 



In principle, to give initial conditions we should specify, at a given time after 
neutrino decoupling, the values of all the multipoles of the neutrino phase 
space distribution. In practice, if the initial conditions are set deep in the 
radiation epoch, the relevant variables only extend, for the purpose of the 
initial conditions, up to the octupole. Specific examples will be given in a 
moment. 

CDM component 

The CDM component is in some sense, the easier. In the standard case the 
evolution equations do not contain neither the magnetic field contribution 
nor the anisotropic stress. The evolution of the density contrast and of the 
peculiar velocity are simply given, in Fourier space, by the following pair of 
equations: 

S' c = tol/-O c , (153) 

e' c +ne c = k 2 <p. (154) 

Magnetized adiabatic and non-adiabatic modes 

The evolution equations of the fluid and metric variables will now be solved 
deep in the radiation-dominated epoch and for wavelengths much larger than 
the Hubble radius, i.e. \kr\ <C 1. In the present lecture only the magnetized 
adiabatic mode will be discussed. However, the treatment can be usefully 
extended to the other non-adiabatic modes. For this purpose we refer the in- 
terested reader to |136] (see also |143j ). Moreover, since this lecture has been 
conducted within the conformally Newtonian gauge, there is no reason to 
change. However, it should be noticed that fully gauge-invariant approaches 
are possible [137] . To give the flavour of the possible simplifications obtainable 
in a gauge-invariant framework we can just use gauge-invariant concepts to 
classify more precisely the adiabatic and non adiabatic modes. For this pur- 
pose, in agreement with Eq. (|126p . let us define the gauge-invariant density 
contrasts on uniform curvature hypersurfaces for the different species of the 
pre-decoupling plasma: 

C 7 = -^+^, C„ = -^ + y, (155) 

Cc = -^ + y, Cb = -^ + y- (156) 

In terms of the variables of Eqs. (|155|) (|156[) the evolution equations for the 
density contrasts, i.e. Eqs. (|144p . (|148[) , (|154[) and (j!54f) . acquire a rather 
symmetric form: 
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From Eqs. (|157|) and (|158[) we can easily deduce a rather important property 
of fluid mixtures: in the long wavelength limit the relative fluctuations in 
the specific entropy are conserved. Consider, for instance, the CDM-radiation 
mode. In this case the non vanishing entropy fluctuations are 

5 7C = -3(C 7 - Cc), = -3(C„ - Cc)- (159) 

Using Eqs. (|157|) and (|158|) the evolution equations for <S 7C and S uc can be 
readily obtained and they are 

s' lc = -{e lh ~e c ), S ' vc = -(e u -9 c ). (ieo) 

Outside the horizon the divergence of the peculiar velocities is C(|fcr| 2 ), so the 
fluctuations in the specific entropy are approximately constant in this limit. 
This conclusion implies that if the fluctuations in the specific entropy are zero, 
they will still vanish at later times. Such a conclusion can be evaded if the 
fluids of the mixture have a relevant energy-momentum exchange or if bulk 
viscous stresses are present [1481 1149] . 

A mode is therefore said to be adiabatic iff £7 = = £ c = Cb- Denoting 
by Q and (j two generic gauge-invariant density contrasts of the fluids of the 
mixture, we say that the initial conditions are non-adiabatic if, at least, we 
can find a pair of fluids for which £j 7^ £j • 

As an example, let us work out the specific form of the magnetized adia- 
batic mode. Let us consider the situation where the Universe is dominated by 
radiation after weak interactions have fallen out of thermal equilibrium but 
before matter-radiation equality. This is the period of time where the initial 
conditions of CMB anisotropics are usually set both in the presence and in the 
absence of a magnetized contribution. Since the scale factor goes, in conformal 
time, as a(r) ~ r and 7i ~ r _1 , Eq. (|9"0|) can be solved for \kr\ <C 1. The 
density contrasts can then be determined, in Fourier space, to lowest order in 
fcr as: 

5-y = S v = —2(f)[ — i? 7 i?B, 

5 h = 6 C = - -RyQu, (161) 

where the fractional contribution of photons to the radiation plasma, i.e. i? 7 
has been introduced and it is related to i.e. the fractional contribution of 
massless neutrinos, as 

r 

1 + r 

r =H^r s °- 68i (f)- <i62) 

In Eq. (|16ip <fii(k) denotes the initial value of the metric fluctuation in Fourier 
space. It is useful to remark that we have treated neutrinos as part of the 
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radiation background. If neutrinos have a mass in the meV range, they are 
nonrelativistic today, but they will be counted as radiation prior to matter- 
radiation equality. Concerning Eq. (|101|) the last remark is that, of course, 
we just kept the lowest order in \kr\ < 1. It is possible, however, to write the 
solution to arbitrary order in \kr as explicitly shown in Ref. |143j . 

Let us then write Eq. (|105|) in Fourier space and let us take into account 
that the background is dominated by radiation. The neutrino quadrupole is 
then determined to be 



R*y k T 



R u 



QR U 



(163) 



where tpi (k) is the initial (Fourier space) value of the metric fluctuation defined 
in Eq. 

Let us then look for the evolution of the divergences of the peculiar veloci- 
ties of the different species. Let us therefore write Eqs. (|146p . (|149|1 and (I153|) 
in Fourier space. By direct integration the following result can be obtained: 



# 7 b 



k t , 



R v f2 B - 4o- B ], 



Rjfij: 



k 2 r-^-a B , 



(164) 
(165) 

(166) 



As a consistency check of the solution, Eqs. (|164|) , (|165|) and (I166P can be 
inserted into Eq. (PTj) . Let us therefore write Eq. (fDT|) in Fourier space 



k 2 H<t>i = AirGa 2 



4 4 

-p 7 (l + p h )dyb + -^pu &u + Pc 



(167) 



where we used that — and we also used the tight-coupling approximation 
since 9 7 = 8h = 0-yb- Notice that in Eq. |9T|) the term arising from the Poynt- 
ing vector has been neglected. This approximation is rather sound within the 
present MHD treatment. In Eq. (|167j) R b < 1 (see Eq. (|147[) for the definition 
of i?b) since we are well before matter-radiation equality. The same observa- 
tion can be made for the CDM contribution which is negligible in comparison 
with the radiative contribution provided by photons and neutrinos. Taking 
into account these two observations we can rewrite Eq. (|167p as 



k 2 H<t>i = 2H 2 {R i e ih + R„e„), 



(168) 



where Eqs. (|97| and |98|) have been used. Inserting then Eqs. (|164j) and (I165|) 
into Eq. (|168p it can be readily obtained that the left hand side exactly equals 
the right hand side, so that the momentum constraint is enforced. 

The final equation to be solved is the one describing the evolution of 
the anisotropic stress, i.e. Eq. (|150p . Inserting Eqs. (|163p and (|165p into Eq. 
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(|150|) we do get an interesting constraint on the initial conditions on the two 
longitudinal fluctuations of the geometry introduced in Eqs. (|89")1 , namely: 

Concerning the magnetized adiabatic mode the following comments are in 
order: 

• the peculiar velocities are always suppressed, with respect to the other 
terms of the solution, by a factor \kr\ which is smaller than 1 when the 
wavelength is larger than the Hubble radius; 

• in the limit ctb — * and — * the magnetized adiabatic mode presented 
here reproduces the well known standard results (see for instance [142] ): 

• the difference between the two longitudinal fluctuations of the metric is 
due, both to the presence of magnetic and fluid anisotropic stresses; 

• the longitudinal fluctuations of the geometry are both constant outside 
the horizon and prior to matter-radiation equality; this result still holds 
in the presence of a magnetized contribution as it is clearly demonstrated 
by the analytic solution presented here. 

The last interesting exercise we can do with the obtained solution is to 
compute the important variables 1Z and £ introduced, respectively, in Eqs. 
(|116|) and l|115jl . Since both if) and <fi are constants for \kr\ < 1 and for 
t < T cq , also 1Z will be constant. In particular, by inserting Eq. (|169|) into Eq. 
(|116|) . the following expression can be obtained: 

Ki = -l(l+ ^R^j <pi - ^(4c7b - Ruthi), (170) 

where TZi(k) denotes the initial value, in Fourier space, of the curvature per- 
turbations. In numerical studies it is sometimes useful to relate the initial 
values of <fr and ip, i.e. <pi and ipi to TZ{. This relation is expressed by the fol- 
lowing pair of formulae that can be derived by inverting Eq. (11701) and by 
using Eq. fi"69"]) : 

10 2i? 7 (4(iB-i? 7 ^B) 



15 + 4i?„ 15 + 4i?„ 

il ~ - 2 i5Tm; ni ~ 5 15 + 4^ {A(JB r ^ b) - (171) 

From the Hamiltonian constraint written in the form (|117p it is easy to deduce, 
in the limit \kr\ <C 1 that d(k) — IZi(k) The same result can be obtained 
through a different, but also instructive, path. Consider the definition of £ 
given either in Eq. (| 1 1 2|) or ()115jl . The variable ( can be expressed in terms of 
the partial density contrasts defined in Eqs. ()155|1 and (|156|) . More precisely, 
from the definitions of the two sets of variables it is easy to show that 
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c 



P'uCv + P'~f(y + PcCc + PbCb 



Pt 



+ Ce 



3(pt + Pt) 



(172) 



Thus, to obtain £ it suffices to find £71 d,, Cb and £c evaluated at the initial 
time and on the adiabatic solution. Using Eqs. (|16ip and p69[) into Eqs. (|155|) 
and (|156p we obtain, as expected, 



C-y = Cu = Cc = Cb = - V>i 



7?. 



1-Sh 



(173) 



This result was expected, since, as previously stressed, for the adiabatic mode 
all the partial density contrasts must be equal. Inserting now Eq. (|173p into 
Eq. (|172[) and recalling that the CDM and baryon contributions vanish deep 
in the radiation epoch, we do get 



c = - 



(174) 



where the last equality follows from the definition of (|116[) evaluated deep in 
the radiation epoch and for the adiabatic solution derived above. 

Up to now, as explained, attention has been given to the magnetized adi- 
abatic mode. There are, however, also other non adiabatic modes that can 
enter the game. We will not go, in this lecture, through the derivation of the 
various non-adiabatic modes. It is however useful to give at least the result 
in the case of the magnetized CDM-radiation mode. In such a case the full 
solution to the same set of equations admitting the adiabatic solutions can be 
written as For the case of the CDM-radiation mode the solution, in the limit 
t < T\ and fcr < 1 can be written as 



ip = ipi\ — 



5y = 

S c = 



5* + — Ryf2ft 



7 7 b 
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k 2 n 
fc 2 n 



fc 2 n 
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R 



"5 I<7 & + -^ L W , i-0i)( — 



where 



Rv 



01 



6i?„ 



15 + AR V 



8(15 
15- 



-2i? w ) 
4R V 



n 



-i? 7 i?B 



(175) 



-i? 7 J?B 

4 ' 



(176) 



r 8(15 + 2i?„) 

In Eq. (|175[) the following notation for the non-vanishing entropy fluctuations 
has been employed: 

S C j = S cv — Sr. {177) 

In deriving Eq. (|175|) it is practical to use a form of the scale factor (ob- 
tained by solving Eqs. (|9~7| , (|9"5|) and for a mixture of matter and radia- 
tion) which explicitly interpolates between a radiation-dominated regime and 
a matter-dominated regime: 



a(r) 



•2(1- 

n 



1 + z oa = 



(178) 



where .f2 m o and J? r o are evaluated at the present time and the scale factor is 
normalized in such a way that ag = 1. In Eq. (|178p 7~i = (2/ Hq) y*' a cq / {2 m o- 
In terms of T\ the equality time is 



r oq = (V2- l)n = H9.07 



/i i7 m o 
0.134 



Mpc, 



(179) 



i.e. 2r cq ~ ri. In this framework the total optical depth from the present to 
the critical recombination epoch, i.e. 800 < z < 1200 can be approximated 
analytically, as discussed in [150] . By defining the redshift of decoupling as 
the one where the total optical depth is of order 1, i.e. >t(zdcc,0) ~ 1, we will 
have, approximately 



1139 



0.0431 / 



Oil 



0.0268 



(180) 



0.6462 + 0.1125 In (J? b /0.0431) ' 

where h = 0.73. From Eqs. (fT80|) and (fT78|) it follows that for 1100 < z dcc < 
1139, 275 Mpc < T doc < 285 Mpc. 

Equations (|178p and (|179p will turn out to be relevant for the effective 
numerical integration of the brightness perturbations which will be discussed 
later on. For numerical purposes the late-time cosmological parameters will 
be fixed, for a spatially flat Universe, as 



w 7 = 2.47 x 10" 



uj h = 0.023, w c = 0.111, u m = u; b +w c , (181) 

where ujx — h 2 S7x and Oa = 1 — ^m] the present value of the Hubble param- 
eter Hq will be fixed, for numerical estimates, to 73 in units of km/ (sec Mpc). 



The values of the cosmological parameters introduced in Eq. (|181[) are compatible 
with the ones estimated from WMAP-3 [1311 11511 1152] in combination with the 
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Transfer matrix and Sachs- Wolfe plateau 

Before presenting some numerical approaches suitable for the analysis of mag- 
netized CMB anisotropics it is useful to discuss a class of analytical estimates 
that allow the calculation of the so-called Sachs- Wolfe plateau. The idea, in 
short, is very simple. We have the evolution equation for £ given in Eq. (|136p . 
This evolution equation can be integrated across the matter-radiation transi- 
tion using the interpolating form of the scale factor proposed in Eq. (|178p . 

Consider, first, the case of the magnetized adiabatic mode where 5p n&< \ = 0. 
Deep in the radiation-dominated epoch, for r -C r cq , — > 1/3 and, from Eq. 
11311), C = 0, so that 

C = Ci - fti, Ci = -\<k (l + j^R^j - ^(^b - R.O B ). (182) 

When the Universe becomes matter-dominated, after r oq , — > and the 
second term at the right hand side of Eq. (I136j) does contribute significantly at 
decoupling (recall that for /i 2 Anatter = 0.134, Td cc = 2.36 r cq ). Consequently, 
from Eq. (|136p . recalling that c 2 = 4a oq /[3(3a + 4a oq )], we obtain 

Cf = Ci- ^ aJ ^ B v Bf = Q Bi . (183) 
4(3a + 4a cq ) 

The inclusion of one (or more) non-adiabatic modes changes the form of Eq. 
(1136|) and, consequently, the related solution (|183p . For instance, in the case 
of the CDM-radiation non-adiabatic mode the relevant terms arising in the 
sum (|133|) are S CJ = <S C „ = Si where Si is the (constant) fluctuation in the 
relative entropy density initially present (i.e. for r <C r oq ). If this is the case 
dpnad — c sPcSi and Eq. (|136|) can be easily solved. The transfer matrix for 
magnetized CMB anisotropies can then be written as 



(184) 
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In the case of a mixture of (magnetized) adiabatic and CDM-radiation modes, 

1 . . 



we find, for a > a oq 



M c< -» 1, M cs -» -g, M CB ■ ( . 
Mss -» 1, Msb -> 0, (185) 



"Gold" sample of SNIa [153] consisting of 157 supernovae (the furthest being 
at redshift z = 1.75). We are aware of the fact that WMAP-3 data alone seem 
to favour a slightly smaller value of u; m (i.e. 0.126). Moreover, WMAP-3 data 
may also have slightly different implications if combined with supernovae of the 
SNLS project |154| . The values given in Eq. (|18ip will just be used for a realistic 
numerical illustration of the methods developed in the present investigation. 
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and A^bb — ¥ 1- Equations (| 1 84|) and (|185p may be used, for instance, to obtain 
the magnetized curvature and entropy fluctuations at photon decoupling in 
terms of the same quantities evaluated for r <C T cq . A full numerical analysis 
of the problem confirms the analytical results summarized by Eqs. (|184p and 
(I185|) . The most general initial condition for CMB anisotropics will then be 
a combination of (correlated) fluctuations receiving contribution from <5p n ad 
and from the fully inhomogeneous magnetic field. To illustrate this point, the 
form of the Sachs- Wolfe (SW) plateau in the sudden decoupling limit will now 
be discussed. 

To compute the SW contribution we need to solve the evolution equation 
of the monopolc of the temperature fluctuations in the tight coupling limit, 
i.e. from Eqs. fUS]) and (fMH]) . 

(186) 

In the sudden decoupling approximation the visibility function, i.e. /C(r) = 
K'(T)e~ K ( T ' and the optical depth, i.e. e _K ( T ' are approximated, respectively, 
by S(t — Tdcc) and by 8(t — Td cc ) (see |155l 1156] for an estimate of the width 
of the last scattering surface). The power spectra of S and f2 B are given, 
respectively, by: 

(, n. n r — 1 / 7 \ n,- 1 

Y J > r s (k)=As\J-j , (187) 

Vn{k)^T{e)Ul L ^-) 2 \ ( 188 ) 
where .4^ , As and Q B l are constants and 



- p ^a\ T ) Pl (T). (189) 

To deduce Eqs. (|187|1 . (|188| and (1189P the magnetic field has been regularized, 
according to a common practice [1281 11251 1126] , over a typical comoving scale 
L = 27r/fci with a Gaussian window function and it has been assumed that 
the magnetic field intensity is stochastically distributed as 

27T 2 

(B l (k 1 T)B3( Pl T)) = —P?(k)P B (k 1 T)5^(k+p), (190) 

where 

?Hk) = (S{ - P B (k,r) = Ab^J- (191) 



B L 



4(6 - e)(2ir) 2s 
e{3-2e)r 2 (e/2)' 
Pbl 



Pb l 



8tt ' 
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As a consequence of Eq. (|190[) the magnetic field does not break the spatial 
isotropy of the background geometry. The quantity k p appearing in Eqs. (|187|) 
and (|19ip is conventional pivot scale that is 0.05 Mpc(see 1132J 11331 1134] for 
a discussion of other possible choices). Equations (|188p and (|189[) hold for 
< s < 1. In this limit the Vo(k) (see Eq. (jl88[) ^ is nearly scale-invariant 
(but slightly blue) . This means that the effect of the magnetic and thermal 
diffusivity scales (related, respectively, to the finite value of the conductivity 
and of the thermal diffusivity coefficient) do not affect the spectrum [126J. 
In the opposite limit, i.e. e 1 the value of the mode-coupling integral 
appearing in the two-point function of the magnetic energy density (and of 
the magnetic anisotropic stress) is dominated by ultra-violet effects related to 
the mentioned diffusivity scales [H>]. Using then Eqs. (fl87 ]l .(fl88 |) and (fT89]) 
the Ce can be computed for the region of the SW plateau (i.e. for multipoles 
£ < 30): 



Ce = 



A 



9 



>5 Zl ( ?v '^ + ym^bl^M) - ^v / AA^ 1 K s ^)cos 7r 



4 3 — — 

-gg AsZ\(n s ,l) - — V Aq R 1 i^BL Z 3 (n r ,e 7 £) cosj br 

_6_ 

25 



As R-t &b LZ 3 (n s ,e, £) cos jbs 



(192) 



where the functions Z\, i? 2 and Z 3 



Z!(nJ) 



7T 2 / k 



4 U, 



n _i r(3-n)r[£+^ 



r>\i-i )rli + |-f 



1 V^/ n(l-e)r(£ + 2-e) 



(193) 



(194) 



e I n 1 

2 ' 4 4 



r2( 1 — £ — » I v I H j- p — £ — a 

1 1 4 2 4lM4^^ 2 4 



(195) 



are defined in terms of the magnetic tilt e and of a generic spectral index 
n which may correspond, depending on the specific contribution, either to 
n r (adiabatic spectral index), or to n s (non-adiabatic spectral index) or even 
to n rs — (n r + n s )/2 (spectral index of the cross-correlation). In Eq. (1192[) 
Irs, Ibr and 7 S {, are the correlation angles. In the absence of magnetic and 
non-adiabatic contributions and for Eqs. (|192j) and Eq. (|193j) imply that for 
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n r = 1 (Harrison-Zeldovich spectrum) £(£ + l)C e /2n = _4 c /25 and WMAP 
data 130J would imply that = 2.65 x 10 -9 . Consider then the physical sit- 
uation where on top of the adiabatic mode there is a magnetic contribution. If 
there is no correlation between the magnetized contribution and the adiabatic 
contribution, i.e. 7^ = tt/2, the SW plateau will be enhanced in comparison 
with the case when magnetic fields are absent. The same situation arises when 
the two components are anti-correlated (i.e. cos7t, r < 0). However, if the fluc- 
tuations are positively correlated (i.e. cos7b r > 0) the cross-correlation adds 
negatively to the sum of the two autocorrelations of C and J?b so that the 
total result may be an overall reduction of the power with respect to the case 
7 br = 7r/2. In Eq. (| 193|) . (| 194[) and (|195[) k = t^ 1 where r is the present 
observation time. 

4.4 Numerical analysis 

The main idea of the numerical analysis is rather simple. Its implementation, 
however, may be rather complicated. In order to capture the simplicity out of 
the possible complications we will proceed as follows. We will first discuss a 
rather naive approach to the integration of CMB anisotropics. Then, building 
up on this example, the results obtainable in the case of magnetized scalar 
modes will be illustrated. 

Simplest toy model 

Let us therefore apply the Occam razor and let us consider the simplest situ- 
ation we can imagine, that is to say the case where 

• magnetic fields are absent; 

• neutrinos are absent; 

• photons and baryons are described within the tight-coupling approxima- 
tion to lowest order (i.e. ctt — * 00); 

• initial conditions are set either from the adiabatic mode or from the CDM- 
radiation mode. 

This is clearly the simplest situation we can envisage. Since neutrinos are 
absent there is no source of anisotropic stress and the two longitudinal fluc- 
tuations of the metric are equal, i.e. <j> — tp. Consequently, the system of 
equations to be solved becomes 



TV 



k 2 c 2 s H 



Pt+Pt 



n 



Spnad, 



(196) 




(197) 



(198) 
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6' c = 3V/ - e c , (200) 

e' c = -w c + fcV- (201) 

We can now use the explicit form of the scale factor discussed in Eq. (|178p 
which implies: 

n _ 1 2(3 + 1) 



T\ x(x + 2) ' 
2 x 2 + 2x 



l x 2 (x + 2) 2 



r 



2 , 1 2(3a; 2 + 6a; + 4) 
W ~ H = ^F x*(x + 2f ' (2 ° 2) 

where x = t/t\. With these specifications the evolution equations given in 
fTgg J) -(|2"DT |) become 



dTZ _ 4x(x + l)(x + 2) i .> 



da; 3 (3x 2 + 6x + 4) 2 



(203) 



_ 3x 2 + 6x + 4 5a; 2 + 10a; + 6 

da 7 ~ ~x(x + l)(x + 2) ~ a; (a; + l)(x + 2)^' ( j 

rf£ 7 _ 4(3x 2 + 6a; + 4) 4(5x 2 + lOx + 6) , 4 - 
dx ~ x(x + l)(x + 2)^ x(x + l)(x + 2) V 3 7b ' 1 j 

<5 7 + k 2 iP, (206) 



dfl 7b _ 2gb (x + 1) « 2 a 

dx i? b + lx(x + 2) 4(1 + i? b ) 7 K 

3(3, 2 + 6a: + 4)^_3(5, 2 + 10, + 6K_^ ^ 



da; x(x + l)(x + 2) x(x+l)(x + 2) 

+ k 2 ^. (208) 



d<9 c 2(x + 1) - 2 



da; x(x + 2) 
In Eqs. (j203[) (j208|) the following rescalings have been used: 

n = kn, e lh = n6 lh , e c = Tl e c . (209) 

The system of equations (|203l) - (|208[) can be readily integrated by giving initial 
conditions for at x\ -C 1. In the case of the adiabatic mode (which is the one 
contemplated by Eqs. (|203[) (|208|) since we set Sp nac i — 0) the initial conditions 
are as follows 

2 

K(xi) = "72.*, i/j(xi) = — -72.*, 
5 7 (xi) = -2^*, 0~ ( b{xi) = 0, 

^c(xi) = -^*, ^c(xi) = 0. (210) 
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It can be shown by direct numerical integration that the system (|203[) - (|208[) 
gives a reasonable semi-quantitative description of the acoustic oscillations. To 
simplify initial conditions even further we can indeed assume a flat Harrison- 
Zeldovich spectrum and set 7?.» = 1. 

The same philosophy used to get to this simplified form can be used to 
integrate the full system. In this case, however, we would miss the impor- 
tant contribution of polarization since, to zeroth order in the tight-coupling 
expansion, the CMB is not polarized. 



Integration of brightness perturbations 

To discuss the polarization, we have to go (at least) to first-order in the 
tight coupling expansion [1571 11581 1159) . For this purpose, it is appropriate 
to introduce the evolution equations of the brightness perturbations of the 
/, Q and U Stokes parameters characterizing the radiation field. Since the 
Stokes parameters Q and U are not invariant under rotations about the axis 
of propagation the degree of polarization P = (Q 2 + U 2 ) 1 / 2 is customarily 
introduced [1591 1160] . The relevant brightness perturbations will then be de- 
noted as A\, Ap. This description, reproduces, to zeroth order in the tight 
coupling expansion, the fluid equations that have been presented before to set 
initial conditions prior to equality. For instance, the photon density contrast 
and the divergence of the photon peculiar velocity are related, respectively, to 
the monopole and to the dipole of the brightness perturbation of the intensity 
field, i.e. <5 7 = 4Z\io and 9 1 = 3k An. The evolution equations of the brightness 
perturbations can then be written, within the conventions set by Eq. (|89|) 



A[ + (ikfi + k')Ai + ikfuj) = tp' + k' 



Aio + nv h - ^P 2 (n)S P 



(211) 



A' P + {ikfi + k')Ap = - [1 - P 2 (fJ,)]S P , (212) 

v' h + Hv h + ikcj)+^- [fie - 4(7 B ] + -£-(«b + 3iAi) = 0. (213) 
4-Kb -Kb 

Equation (|213p is nothing but the second relation obtained in Eq. (|140|) having 
introduced the quantity ikv\> = 9^. The source terms appearing in Eqs. (|21ip 
and (|212p include a dependence on P2{n) = (3^ 2 — 1)/2 ( Pi{n) denotes, in this 
framework, the £-th Legendre polynomial); ; fi — k ■ n is simply the projection 
of the Fourier wave-number on the direction of the photon momentum. In 
Eqs. (|21ip and (|212p the source term Sp is defined as 

S P {k, r) = A I2 (k, r) + Ap (k, r) + A P2 (k, r). (214) 

The evolution equations in the tight coupling approximation will now be 
integrated numerically. More details on the tight coupling expansion in the 
presence of a magnetized contribution can be found in 136\. 
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The normalization of the numerical calculation is enforced by evaluating, 
analytically, the Sachs- Wolfe plateau and by deducing, for a given set of spec- 
tral indices of curvature and entropy perturbations, the amplitude of the power 
spectra at the pivot scale. Here is an example of this strategy. The Sachs- Wolfe 
(SW) plateau can be estimated analytically from the evolution equation of 1Z 
(or £) by using the technique of the transfer matrix appropriately generalized 
to the case where, on top of the adiabatic and non-adiabatic contributions 
the magnetic fields are consistently taken into account. The main result is 
expressed by Eq. (|192[) . 

If the SW plateau is determined by an adiabatic component supplemented 
by a (subleading) non-adiabatic contribution both correlated with the mag- 
netic field intensity the obtainable bound may not be so constraining (even 
well above the nG range) due to the proliferation of parameters. A possible 
strategy is therefore to fix the parameters of the adiabatic mode to the values 
determined by WMAP-3 and then explore the effect of a magnetized contri- 
bution which is not correlated with the adiabatic mode. This implies, in Eq. 
(|192j) that As = and 7b r = tt/2. Under this assumption, in Figs. and 
[9] the bounds on are illustrated. The nature of the constraint depends, 
in this case, both on the amplitude of the protogalactic field (at the present 
epoch and smoothed over a typical comoving scale L = 2-k jkj}) and upon its 
spectral slope, i.e. e. In the case e < 0.5 the magnetic energy spectrum is 
nearly scale-invariant. In this case, diffusivity effects are negligible (see, for 
instance, [T91I125] ). As already discussed, if e ^> 1 the diffusivity effects (both 
thermal and magnetic) dominate the mode-coupling integral that lead to the 
magnetic energy spectrum Q]|l I125j . 

In Fig. [8] the magnetic field intensity should be below the different curves 
if the adiabatic contribution dominates the SW plateau. Different choices of 
the pivot scale k p and of the smoothing scale Ul, are also illustrated. In Fig. [8] 




Fig. 8. Bounds on the protogalactic field intensity as a function of the magnetic 
spectral index e for different values of the parameters defining the adiabatic contri- 
bution to the SW plateau. 
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the scalar spectral index is fixed to n r = 0.951 |131j . In Fig. [S] the two curves 
corresponding, respectively, to n r = 0.8 and n r = 1 are reported. 




-8.6 - 



Fig. 9. Same plot as in Fig. [5] but with emphasis on the variation of n r . 

If e < 0.2 the bounds are comparatively less restrictive than in the case e ~ 
0.9. The cause of this occurrence is that we are here just looking at the largest 
wavelengths of the problem. As it will become clear in a moment, intermediate 
scales will be more sensitive to the presence of fully inhomogeneous magnetic 
fields. 

According to Figs. [5] and [5] for a given value of the magnetic spectral index 
and of the scalar spectral index the amplitude of the magnetic field has to 
be sufficiently small not to affect the dominant adiabatic nature of the SW 
plateau. Therefore Figs.[5]and[n](as well as other similar plots) can be used to 
normalize the numerical calculations for the power spectra of the brightness 
perturbations, i.e. 

^|A(fc,r)| 2 , ^\A P (k,r)\ 2 , ^\Mk,r)A P (k,r)\. (215) 

Let us then assume, for consistency with the cases reported in Figs. [5] and 
[5J that we are dealing with the situation where the magnetic field is not 
correlated with the adiabatic mode. It is then possible to choose a definite 
value of the magnetic spectral index (for instance e = 0.1) and a definite 
value of the adiabatic spectral index, i.e. n r (for instance n r — 0.951, in 
agreement with |131j ). By using the SW plateau the normalization can be 
chosen in such a way the the adiabatic mode dominates over the magnetic 
contribution. In the mentioned case, Fig. [8] implies Bl < 1.14 x 10~ 8 G for a 
pivot scale k p = 0.002 Mpc -1 . Since the relative weight of the power spectra 
given in Eqs. (|187|l and (|188| is fixed, it is now possible to set initial conditions 
for the adiabatic mode according to Eqs. (fT6T|) - (fT63| . pM | - (fl66|l and {Tog]) 
deep in the radiation-dominated phase. The initial time of integration will be 
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chosen as T\ = 10~ 6 Ti in the notations discussed in Eq. (|178p . According to 
Eq. (|179p . this choice implies that t\ <C r cq . 

The power spectra of the brightness perturbations, i.e. Eq. (|215p . can be 
then computed by numerical integration. Clearly the calculation will depend 
upon the values of w m , cub, oj c and R v . We will simply fix these parameters to 
their fiducial values reported in Eqs. (|181| (see also (|147ll ) and we will take 
N v = 3 in Eq. (|162| determining, in this way the fractional contribution of 
the neutrinos to the radiation plasma. 

The first interesting exercise, for the present purposes, is reported in Fig. 
[TOl where the power spectra of the brightness perturbations are illustrated for 
a wave-number k = 0.1 Mpc -1 . Concerning the results reported in Fig. [10] 
different comments are in order: 

• for £ = 0.1 and n r = 0.951, the SW plateau imposes B L < 1.14 x 10~ 8 G; 
from Fig. [TO] it follows that a magnetic field of only 30 nG (i.e. marginally 
incompatible with the SW bound) has a large effect on the brightness 
perturbations as it can be argued by comparing, in Fig. I10| the dashed 
curves (corresponding to 30 nG ) to the full curves which illustrate the 
case of vanishing magnetic fiels; 

• the situation where Bl > nG cannot be simply summarized by saying that 
the amplitudes of the power spectra get larger since there is a combined 
effect which both increases the amplitudes and shifts slightly the phases 
of the oscillations; 

• from the qualitative point of view, it is still true that the intensity oscillates 
as a cosine, the polarization as a sine; 

• the phases of the cross-correlations are, comparatively, the most affected 
by the presence of the magnetic field. 

The features arising in Fig. [TO] can be easily illustrated for other values of e 
and for different choices of the pivot or smoothing scales. The general lesson 
that can be drawn is that the constraint derived only by looking at the SW 
plateau are only a necessary condition on the strength of the magnetic field. 
They are, however, not sufficient to exclude observable effects at smaller scales. 
This aspect is illustrated in the plot at the left in Fig. [ITJ which captures a 
detail of the cross-correlation. The case when Bl = can be still distinguished 
from the case Bl = 0.5 nG. Therefore, recalling that for the same choice of 
parameters the SW plateau implied that Bl < 11.4 nG, it is apparent that 
the intermediate scales lead to more stringent conditions even for nearly scale- 
invariant spectra of magnetic energy density. For the range of parameters of 
Fig. [H] we will have that Bl < 0.5 nG which is more stringent than the 
condition deduced from the SW plateau by, roughly, one order of magnitude. 

If e increases to higher values (but always with e < 0.5) by keeping fixed 
Bl (i.e. the strength of the magnetic field smoothed over a typical length 
scale L = the amplitude of the brightness perturbations gets larger 

in comparison with the case when the magnetic field is absent. This aspect 
is illustrated in the bottom plot of Fig. [Tl] where the logarithm (to base 10) 
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Fig. 10. The power spectra of the brightness perturbations for a typical wave- 
number k = O.lMpc -1 . The values of the parameters are specified in the legends. 
The pivot scale is k p = 0.002 Mpc -1 and the smoothing scale is kL = Mpc -1 (see 
Figs. [8] and O. 




1 0.i 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

e 

Fig. 11. A detail of the cross-correlation (top). The autocorrelation of the intensity 
at T c icc as a function of e, i.e. the magnetic spectral index (bottom). 



of the intensity autocorrelation is evaluated at a fixed, wave-number (and at 
Tdcc) as a function of e. The full line (corresponding to a Bj, = 10 nG) is 
progressively divergent from the dashed line (corresponding to Bl = 0) as e 
increases. 

In Fig. [T2] the power spectra of the brightness perturbations are reported 
at Tdec and as a function of k. In the two plots at the top the autocorrelation 
of the intensity is reported for different values of Bl (left plot) and for dif- 
ferent values of £ at fixed Bl (right plot). In the two plots at the bottom the 
polarization power spectra are reported always at Td cc and for different values 
of Bl at fixed e. The position of the first peak of the autocorrelation of the 
intensity is, approximately, fed — 0.017 Mpc -1 . The position of the first peak 
of the cross-correlation is, approximately, 3/4 of fed- From this consideration, 
again, we can obtain that Bl < 0.3 nG which is more constraining than the 
SW condition. 

Up to now the adiabatic mode has been considered in detail. We could 
easily add, however, non-adiabatic modes that are be partially correlated with 
the adiabatic mode. It is rather plausible, in this situation, that by adding 
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new parameters, also the allowed value of the magnetic field may increase. 
Similar results can be achieved by deviating from the assumption that the 
magnetic field and the curvature perturbations are uncorrelated. This aspect 
can be understood already from the analytical form of the SW plateau (|192p . If 
there is no correlation between the magnetized contribution and the adiabatic 
contribution, i.e. 7& r = tt/2, the SW plateau will be enhanced in comparison 
with the case when magnetic fields are absent. The same situation arises 
when the two components are anti-correlated (i.e. cos7& r < 0). However, if 
the fluctuations are positively correlated (i.e. cos7& r > 0) the cross-correlation 
adds negatively to the sum of the two autocorrelations of TZ and J?b so that 
the total result may be an overall reduction of the power with respect to the 
case jbr = n/2. 



k p ^ 0.002 Mpc"' k^l Mpc~ ! ,n r *0.951.E^0.1 ^ k^ = 0.002 Mpc" 1 . k^ 1 Mpc"' ri f ^ 0.951, B L = 10 riG 




log{k Mpc) log{k Mpc) 



Fig. 12. The power spectra of the brightness perturbations at Td cc for the parame- 
ters reported in the legends. 

From Fig. [T2l various features can be appreciated. The presence of magnetic 
fields, as already pointed out, does not affect only the amplitude but also the 
phases of oscillations of the various brightness perturbations. Moreover, an 
increase in the spectral index e also implies a quantitative difference in the 
intensity autocorrelation. 
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5 Concluding remarks 

There is little doubts that large-scale magnetic exist in nature. These fields 
have been observed in a number of different astrophysical systems. The main 
question concerns therefore their origin. String cosmological models of pre-big 
bang type still represent a viable and well motivated theoretical option. 

Simple logic dictates that if the origin of the large-scale magnetic fields 
is primordial (as opposed to astrophysical) it is plausible to expect the pres- 
ence of magnetic fields in the primeval plasma also before the decoupling 
of radiation from matter. CMB anisotropies are germane to several aspect 
of large-scale magnetization. CMB physics may be the tool that will finally 
enable us either to confirm or to rule out the primordial nature of galactic 
and clusters magnetic fields seeds. In the next five to ten years the forth- 
coming CMB precision polarization experiments will be sensitive in, various 
frequency channels between 30 GHz and, roughly 900 GHz. The observations 
will be conducted both via satellites (like the Planck satellite) and via ground 
based detectors (like in the case of the QUIET arrays). In a complementary 
view, the SKA telescope will provide full-sky surveys of Faraday rotation that 
may even get close to 20 GHz. 

In an optimistic perspective the forthcoming experimental data together 
with the steady progress in the understanding of the dynamo theory will hope- 
fully explain the rationale for the ubiquitous nature of large-scale magnetiza- 
tion. In a pessimistic perspective, the primordial nature of magnetic seeds will 
neither be confirmed nor ruled out. It is wise to adopt a model-independent 
approach by sharpening those theoretical tools that may allow, in the near 
future, a direct observational test of the effects of large-scale magnetic fields 
on CMB anisotropies. Some efforts along this perspective have been reported 
in the present lecture. In particular, the following results have been achieved: 

• scalar CMB anisotropies have been described in the presence of a fully 
inhomogeneous magnetic field; 

• the employed formalism allows the extension of the usual CMB initial 
conditions to the case when large-scale magnetic fields are present in the 
game; 

• by going to higher order in the tight coupling expansion the evolution of 
the brightness perturbations has been computed numerically; 

• it has been shown that the magnetic fields may affect not only the ampli- 
tude but also the relative phases of the Doppler oscillations; 

• from the analysis of the cross-correlation power spectra it is possible to 
distinguish, numerically, the effects of a magnetic field as small as 0.5 nG. 

It is interesting to notice that a magnetic field in the range 10 _10 -10~ n G 
is still viable according to the present considerations. It is, therefore, not ex- 
cluded that large-scale magnetic fields may come from a primordial field of the 
order of 0.1-0.01 nG present prior to gravitational collapse of the protogalaxy. 
Such a field, depending upon the details of the gravitational collapse may be 
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amplified to the observable level by compressional amplification. The present 
problems in achieving a large dynamo amplification may therefore be less rel- 
evant than for the case when the seed field is in the range 10 _9 nG = 10~ 18 nG. 
To confirm this type of scenario it will be absolutely essential to introduce the 
magnetic field background into the current strategies of parameter extraction. 

The considerations reported in the present lecture provide already the 
framework for such an introduction. In particular, along a minimalist per- 
spective, the inclusion of the magnetic field background boils down to add 
two new extra-parameters: the spectral slope and amplitude of the magnetic 
field (conventionally smoothed over a typical comoving scale of Mpc size) . The 
magnetic field contribution will then slightly modify the adiabatic paradigm 
by introducing, already at the level of initial conditions, a subleading non- 
Gaussian (and quasi-adiabatic) correction. 
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